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It is shown that the effective interaction strength of three bosons at small collision energies can be extracted 
from their wave function at zero energy. An asymptotic expansion of this wave function at large interparticle 
distances is derived, from which is defined a quantity D named three-body scattering hypervolume, which is 
an analog of the two-body scattering length. Given any finite-range interaction potentials, one can thus predict 
the effective three-body force from a numerical solution of the Schrodinger equation. In this way the constant 
D for hard-sphere bosons is computed, leading to the complete result for the ground state energy per particle 
of a dilute Bose-Einstein condensate (BEC) of hard spheres to order p 2 , where p is the number density. Effects 
of D are also demonstrated in the three-body energy in a finite box of size L, which is expanded to the order 
L~ 7 , and in the three-body scattering amplitude in vacuum. Another key prediction is that there is a violation of 
the effective field theory (EFT) in the condensate fraction in dilute BECs, caused by short-range physics. EFT 
predictions for the ground state energy and few-body scattering amplitudes, however, are corroborated. 



I. INTRODUCTION 

The ground state energy per particle of a dilute Bose 
Einstein condensate (BEC) is 
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plus higher order terms in number density p, where a is the 
scattering length, w = Att/3— \/3 — 2.4567 • • • , and £ 3 is a 
constant. The first three terms in this expansion were discov- 
ered, respectively, by Bogoliubov [1], Lee, Huang, and Yang 
[2, 3], and Wu and others [4-6]. 

The £3 term has remained the least understood. It was 
known to Wu [4] that £3 is given by a parameter £3 for the 
ground state energy of 3 bosons in a periodic cubic volume 
[4], plus many-body corrections. Braaten and Nieto fully de- 
termined these many-body corrections using the effective field 
theory (EFT) [7] but, like Wu, they left undetermined a param- 
eter <?3 (k) for the effective strength of the three-body force 
near the scattering threshold [7], which is related to Wu's pa- 
rameter £3 [4]. [The difference, g^n) — g^(n'), is known for 
any momentum scales k and k' [7].] 

For bosons with large scattering length, g 3 (k) was recently 
computed using the EFT [8]. For other model interactions, 
Braaten et al [9] used the Monte Carlo (MC) results of the 
energy density [10] to extract £3 but, because of statistical un- 
certainties of the MC data that are difficult to reduce, they did 
not obtain satisfying answers [9]. In summary, g$ (k) remains 
unknown for almost all bosonic systems. 

In this paper this three-body force and a few related proper- 
ties of the TV-body system are studied. Implications for dilute 
Bose-Einstein condensates are also explored. 

We know that the interaction of two bosons at low energy 
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interaction) is dominated by the two-body scattering length 
a, while a is present in the small-momentum expansion of 
the two-boson wave function at zero collision energy and zero 
orbital angular momentum: 

q = (2n) 3 S{q) - Ava/q 2 + u + 0{q 2 ). (2) 
Analogously, the effective interaction strength of three identi- 
cal bosons at low energy should be present in the small-q ex- 



pansion of the wave function of the same three bosons at zero 
collision energy and zero orbital angular momentum, 0j c 3 1 ' k2k3 ■ 
It is shown in this paper that this is indeed the case. For small 
momenta hqi ~ hq (X}i=i 1i = 0), it is found that 

^Lm = (27r) 6 5( qi )5(q 2 ) + G qiq2q3 
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+ 0(q- 1 ), (3) 

= 2j{q\ + q 2 + qf), and that D, named three- 
body scattering hypervolume (with dimension [length] 4 ), is a 
suitable parameter for the three-boson interactional 1] 

The definition of the three-body parameter in Eq. (3) per- 
mits one to determine this parameter in an elementary way, 
namely by solving the three-body Schrodinger equation in 
free space and matching the solution to Eq. (3) at small mo- 
menta or its Fourier transform at large relative distances. 

(3) 

In Sec. II of the present paper, </> qi q 2 q 3 is expanded to the 
order q 1 at small momenta [Eqs. (40)] and, correspondingly, 
its Fourier transform expanded to the order R~ 7 at large rel- 
ative distances R [Eqs. (45)]. Although the parameter D 
first appears in the latter expansion at the order i? -4 , three 
higher order corrections are determined, to facilitate much 
more accurate determinations of D from numerical solutions 
to Schrodinger equation at not-so-large relative distances. 

In Sec. Ill, the results of Sec. II are applied to bosons inter- 
acting through the hard-sphere (HS) potential. By solving the 
3-body Schrodinger equation numerically the author found 



£>hs = (1761.5430 ±0.0024)a 4 
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In Sec. IV the ground state of three identical bosons in a 
large periodic cubic volume of side L is determined perturba- 
tively in powers of L _1 . The energy is expanded [12] to the 
order L~ 7 (h = m boson = 1): 
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where r s is the two-body effective range. 

In Sec. IV B Wu's parameter £ 3 [4] is expressed in terms of 
D; £ 3 for hard-sphere bosons is then found [Eq. (109)]. 

In the rest of Sec. IV the A/"-boson energy and momentum 
distribution are found. The importance of the parameter w 
[defined in Eq. (2)], which is absent in the effective field the- 
ory, is stressed. General formulas for the BEC energy and 
condensate depletion are obtained [Eqs. (118) and (119)]. 

In Sec. V the scattering amplitude of three bosons at low 
energy is computed [Eqs. (135)]. After a discrepancy be- 
tween Ref. [7] and our result at r s ^ is resolved, g 3 (re) 
of Ref. [7] is expressed in terms of the scattering hypervol- 
ume D [Eq. (136)]. The ground state energy per particle of a 
dilute Bose gas of hard spheres is finally determined to order 
p 2 [Eq. (139)]. 

II. ASYMPTOTICS OF (3) AT SMALL MOMENTA OR 
LARGE RELATIVE DISTANCES 

We consider identical bosons with instantaneous interac- 
tions that are translationally, rotationally, and Galilean invari- 
ant, and finite-ranged (ie, limited within a finite interparti- 
cle distance r e ). So the 2-body interaction |C/kik 2 k 3 k 4 con- 
serves momentum, is invariant under rotation or any equal 
shift of k/s (1 < i < 4), and is smooth. Also, because of 
Bose statistics, we can symmetrize U with respect to the in- 
coming (outgoing) momenta without losing generality [13]: 

£/k 1 k 2 k 3 k 4 = ^k 2 k 1 k 3 k 4 = ^k 1 k 2 k 4 k 3 = ^k 4 k 3 k2k! ( the l£lSt 

equality is the Hermiticity condition). Similar properties are 
held by the 3-body interaction, ^U^^^^. 

We will consider scattering states only. Units such that H = 
TO boson = 1 are used. 

A. Two-body special functions 

This subsection is in preparation for the rest of the paper. 
We define 2-body special functions in the /-wave channel 
(I is even for identical bosons), with zero magnetic quantum 



number along the direction n: <^A£, g^, [14] 

(M°)k = 0, (Hf^) k = 4l {HgW) k = fW, (6) 

(HX) k = k 2 X k + l -J ^ Ukk ,X k ,, (7) 

where J7kk' = ^k,-k,k'.-k'- For / = 0, we write these func- 
tions simply as </>k, /k, ,9k, For I > 2, we use usual 

symbols d,g,i, . . . to represent I = 2, 4, 6 Let 

and similarly for the other functions. The amplitudes of these 
functions are fixed by the following equations at r > r e : 
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where Pi is the Legendre polynomial [P/(l) = 1], and the 
Z-wave scattering phase shift Si at low energy k 2 satisfies: 

k 2l+1 cot Si(k) = -a^ +r l k 2 /2\ + r' l k i /4\ + 0{k 6 ), (9) 



where a — a, r — r s , r' = r' s , a 2 = a d , r 2 = r d , a 4 = a g , 
and a 6 = <Zj. Now define harmonic polynomials 

Q ( l\k) = k l Pi(h-k). (10) 

The following more general formulas in the momentum space 
are directly related to Eqs. (8) (by Fourier transformation): 
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where i = \/— 1 7^ i, and Z/fc 4 and Z//c 6 are generalized 
functions (in this paper Z is merely a symbol and not a num- 
ber): 
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derive the following asymptotic expansions: 
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Z/fc 4 and Z/fc 6 have —00 values at k = to cancel certain 
integrals shown above. They inevitably arise from the Fourier 



where 7q* q2q3 and S^'^ both scale like q s (including possi- 
bly q s In" q, n = 1, 2, . . . )• The minimum values of s in the 



transformation of functions such as |r|. The Z-functions are above equations will be justified below, 
more completely described in Appendix A. 



The infinite series such as J2i=o ii^pk 21 in Eqs. (11) ac- 
count for the deviations of <^(r), f n l \r), and g n L \r) from 



Eqs. (8) at r < r e . Since r e < 00, these series are convergent. 

(I) ~(i) 

The superscripts of u\ and f- will be omitted at I = 0. 
From Eqs. (6) and (1 1) we derive that at small k 
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For any unknown X k we have the uniqueness theorem [15]: 
(HX) k = (all k) and X k = o(fc~ 3 ) (small k) => X k = 0. 



The following identity is needed in the analysis of the mo- 
mentum distribution of N particles at low density (K stands 
for the real part) [16]: 
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B. Asymptotics of k 3 lk2k at small momenta 

From this point on, we let q's be small momenta and qi's 
scale like q 1 , while k's be independent from q's. We will 



The 3-body Schrodinger equation can be written in two spe- 
cial forms [17]: 
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J k , is the shorthand for J , and / k , k „ stands for 



(13) JJ^'^" 



At small q's we have Taylor expansions: 
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qiq2q3 

n= E K^< 0) - 

s=0,2,4... 



= K0 + Ki(9 2 +<Z 2 2 +^) + O(g 4 ), (22a) 
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Now we fix the overall amplitude of (f>^ (X^=i 1* = 0) : 
^ q 3 iq 2 q3 = (27r) 6 <5(qi)<5(q2) + higher order terms. Because 
<5(Aq) = \- 3 S(q), <5(q) scales like q- 3 . So s > -6 for T«, 
and 



^3 = (2^) 6 <5(qi)<5(q 2 ). 
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indicating that s > —3 for 

The following statement is now true at si = —6: 
Statement s±: all the functions for s < si, and all the 

for s < si + 2, have been formally determined. 
We can then do the following expansions at small q, for 

— 6 < s < Si. 
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where f q k scales like ^fc" ". Note also that S k " 
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Therefore, the asymptotic expansion 



of S' k si+3 ^ q at small k has been determined to the order k 3 : 

-3 

s (-i+3)q = J- t { X 3 ' m) +0(q Sl+3 k- 2 ). (25a) 

m— — si —9 

Equations (24) and (25 a) ensure the continuity of the wave 
function <j)^ across two connected regions. Extracting all the 
terms that scale like q si + 3 from Eq. (18), we get 

(jff5 ( Sl +3)q )k = _3g2 5 ( S i+l)q /4 _ W £+V q Sl+\ (25b) 

If wj£ +3) ^ (ie, if si + 3 = 0, 2, 4 • • • ), we take it as 
formal input. Solving Eqs. (25), with the help of (6) and (11), 

we thus determine 5 k Sl+3 - )q . The uniqueness of the solution 
is guaranteed by (13). 

Once 5^ Sl+3 ' is determined, the right hand side of Eq. (17) 
can be determined up to the order q Sl+3 , if the coefficients of 
the Taylor expansion of ^q iq2q3 [see Eq. (22a)] are regarded 
as formal input. Solving Eq. (17), and noting that G~* q2q3 <~ 
q 2 , we thus determine T^ Sl+1 \ 

So now the truth of Statement (s\ + 1) is established. 

We can thus formally determine all the functions and 
S^ s \ by repeating the above routine, starting from si = —6. 
The results of this program are shown below. 

Step 1. 4~ 3)q = (27r) 3 <5(q)(2^) 3 (5(k) + 0(<r 3 fc- 2 ) at 
small k, and (HS ( - 3 ^) k = 0, so 



k- 



(26) 



Step 2. With the help of Eqs. (12), and noting that 

Ei=i = °. we s et 



T ^L = E -(^a/p^fSidi). (27) 



( — 5) 

Expanding T q _ q / 2 +k - q /2-k at sma U 1 we 8 et a ^ tne func- 
tions tq k '"^ for n + m = —5; shown below are those in the 
range — 2 < n < 3 (zeros are omitted): 



iqk 2 '~ 3) = -(8»/9 2 )(^) 3 *), 
i q °k" 5) =- 7 ra(q-V fe ) 2 (27r) 3 5(k), 



(28a) 
(28b) 



*qk~ 7) = -(W48)(q • V k ) 4 (27r) 3 <5(k)^. (28c) 



Step 3. s£ 2)q = t qk 2 ' 3) + 0(q- 2 k- 2 ) at small k, and 

(ffS(-2)q) k = o, so 

$(-2)* = _^na/q 2 )<Pu. (29) 
Step 4. Using the same method as in Step 2, we get 



rp{-i) _ or. 2 2^. V^„-2 
J qiq2qa — OZ7r U ^919293 ^ </i ■ 



(30) 



»=1 



At small q we have the following expansions (named "Z-5 
expansions"; see Appendix B for details): 

(k 2 + 'iq 2 /^- 1 = k- 2 - V3(27r) 3 <S(k)g/87r - 3q 2 Z/4k 4 
+ V3V 2 k (2ir) 3 5(k)q 3 /64ir + 9q 4 Z/16k 6 
- 3V3Vt(27r) 3 (5(k)g 5 /51207r + 0(<? 6 ), (31a) 

(|k + q/2|- 2 + |k - q /2|- 2 )/(fc 2 +3 9 2 /4) = (27rf5(k)/6q 
+ 2Z/k 4 + [ - V 2 J 48 + (1/24 - V3/16n)Q^{W k )] 



x (2ir) 3 6(k)q + [4ZQ^(k)/3fc 8 - 4Z/3k 6 ]q 2 
+ [V^/1280 - (1/448 - 3V3/8967r)V^g q d) (V fc ) 
+ (19/10080- 3V3/8967r)g q 9) (V fe )](27r) 3 (5(k)g 3 

+ 0(g 4 ), (31b) 



soT 



(-4) 



k = Er=-2*qT"" ) ' Where 



q,-q/2+k,-q/2- 



4k 2 '~ 2) = ^ 2 a 2 /q 2 k 2 , 
,(-1-3) 

q,k 



(0,-4) 
q,k 



47T«;a 2 (27r) 3 ( 5(k)/( ? , 
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(32a) 
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and for brevity higher order terms [which can readily be ob- 
tained from Eqs. (31)] are not shown. 

Step 5. S* k ~ 1)q = t q ; k '~ 3) + Oiq^k' 2 ) at small k, and 

(#S(-i)q) k = o, so 



= (^wa 2 /q)fa. 
Step 6. This is similar to Steps 2 and 4. 



(33) 



T £ZL = -l^ 2 wa 3 G qiq2q3 + «o ^(27r) 3 <5( qi ). 
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Doing the Z-5 expansion (Appendix B), one finds that 
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(-3) 



q,- q /2+k,-q/2-k = En=-3 *q"k 3 n) ■ For brevity only 
one of the terms is shown here (to be used in the next step): 



t 



(0,-3) 
q,k 



[I6wa 3 ln(qr| |) + 2u + (Un/VS - 16)wa 3 



x (27r) 3 (5(k) - 32Tr 2 wa 3 Z 1/lal (k)/k :i . (35) 



This equation remains unaffected if both \a\'s in the logarithm 
and in the subscript of Z are replaced by any other length scale 
simultaneously, because of Eq. (A10). 
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Step 7. Substituting s\ = —3 into Eq.(25b), we get 
(HS^) k = 67ra0 k - . Because the right hand side 
does not depend on k or q, we can introduce a single function 
ci k , which is independent of k and satisfies 

(Hd) k = 67ra0 k - W^ 0) , (36a) 

and get S^* 1 = <i k + (linear combination of qk ). Equa- 
tion (36a) does not completely determine d k , since d k + r?(/> k 
satisfies the same equation. At small k [Eq. (25a)], 

4° )q = [-2 7 raOf(V fe )/3+16 U ; a 3 ln( (Z | a |)](2^) 3 (5(k) 

+ d a (k)+0(k- 2 ), 

where 

d a (k) = -(7ra/3)V^(27r) 3 <5(k) + 40n 2 a 2 Z / k 4 

+ [2m + (14tt/V3 - 16)wa 3 ] (27r) 3 <5(k) 

- 3277 2 wa 3 Z 1/H (k)/k 3 (36b) 

does not depend on k or q. Noting Eq. (11a), we can now 



complete our defintion of d k [18] and determine . 

d k = d„(k) + 0(k- 2 ) at small k. (36c) 

4° )q = 16wa 3 ln(q\a\)ci> k + KWjk + d k- (37) 

Step 8. This is similar to Steps 2, 4, and 6. Extracting all 
the terms that scale like q° from both sides of Eq. (17), and 
solving the resultant equation, we get 

^qiqs = G ?i92 93 [ - 647r W a 4 m( gig2 <? 3 H 3 ) - D] 

3 

-Y,8Trau /q 2 , (38) 

i=l 

(39) 

The subsequent steps are similar to the above ones but con- 
siderably lengthier; details will not be shown. At the comple- 
tion of Step 14 the following results are accumulated: 
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-&<f ln^laD-Ca? 3 Hq\a\)+^q 



, . Sua 4nwa 2 , , , /- 4 32v / 3wa 5 9, ? 
(27r) 3 <5(q) - — + — - — + 16wa 3 \n(q\a\) + 24:V3wa 4 q \n(q\a\) + &q 

-3nwa 2 q - I2wa 3 q 2 ln(g|a|) - 18V3u;aV ln(g|a|) - ^-q 3 
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g 2 ln 2 ( g |a|) 

'i 

9ttw 9 3 
Jk + — 7-a q g k 



IOtto - 107r(27r - 3\/3)a 2 g - 4ioaV ln(g|a|) + ^V] ^ + ^"( 27r ~ 3\/3)aV4k 
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9tt 



3tt 



^-ag 2 + ^(76tt - 135V3>V 0$ + rf k + g 2 ^ + 0(<z 4 ), (40a) 



dsqs = (2-) 3 ^(qi)(2 7 r) 3 ( 5(q 2 ) + ^ , 



]T [-47ra(27r) 3 5(<fc) 



327r 2 a 2 167r 2 wa 3 4l . , ,, D 

2 bA-Kwa ln(gj|a|) — — 

Qi Qi 3 



96\/3Tvwa 5 qi \n( qi \a\) - Ana^iqi + l28V3wa 6 q 2 ln 2 (^|a|) + Aira&ql Mu\ a \) + 407r 2 aa d Q^ ) (p i ) + Ana^q 3 ln(&|a|) 

3 8 4^ 

- Ana^q 3 - 40tt 2 (2tt - 3 V3) a 2 a d tnQ { ^ (p<)l + ^ [ ( 27r ) 3(5 (^ ) ( u o + "iP 2 + u-jp?) - (i»o + «iP?) + («o + «ip 2 ) 



i=l 



+ 16wa 3 ln(^|a|)«o + 24:V3wa 4 qi \n(qi\a\)u + Ci^uo - 3mva 2 qif + Xo + 0(g 2 ), (40b) 

I 



where q 2 c?S is a quadratic polynomial of q [and for any rota- £3^ = 5(9\/3 — 47r)wa 4 /3 + 15-7r(27r — 3\/3)a 2 adrrf/4, 



*qk 

tion r, cg rk = 4k]' and 



(41c) 



6 = V3D/8ir - 8(V3 - 7r/3)u;a 4 - 37rwa 3 r s /2, (41a) 



6 =aZ?/V3 7T 2 -(260/9+128/\/3 7r)wa 5 +2wa 4 r ;s , (41b) 



C 3 = (24/tt + l6/VS)wa e + 9V3wa 5 r s , (4 Id) 
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£ 3 = - [(1/8tt 2 + l/12V3 7r)a + 3\/3r s /647r]aL> 
+ (17V3/2 + 22/tt + 353tt/27) too 
+ (llVS/4 - 7ir/6)wa 5 r s + 9irwa 4 r 2 jm 
+ 3irwa 3 r' s /32 + (IOtt 2 - 45V3tt/4) aa d , (41e) 

Xo = 97ra Ul - 3wi/2 - 2/ti - y C/ k'4k" ( 42 ) 

ki is defined in Eq. (22a), J k , = / j^jp , and wi is a coeffi- 
cient in the following Taylor expansion at small k: 

^ ^ £/kk'dk' - + ^k 2 + 0(k 4 ). (43) 

The other symbols in Eqs. (40) are defined previously: a, r s , 
r' s , a d , and r d in Eq. (9), </> k , /k, 5k, 4k' and <^qk in 
Eqs. (6) and (11), tt; and /o in Eqs. (11) with Z = 0, w = 
4tt/3 - V3, Eti Qi = 0, p, in Eq. (19), in Eq. (10) 
with I = 2, and d k in Eqs. (36). 

At the order q 2 in the expansion of 0q 3 iq2 q 3 , one 
will encounter another 3-body parameter, D', in the term 
—D'G qiq2 q 3 X)i=i [Another contribution at the same or- 
der, Xi (q 2 + q 2 + q 2 ), is smooth and contributes nothing to the 
Fourier-transformed wave function at large relative distances.] 
D' is independent from D for general interactions. 



In general, if one expands </> qi q 2 q 3 to the order q m , one will 
encounter a total of exactly 

N 3 (m) = round(m' 2 /48 + m'/3 + 89/72) (44) 

independent 3 -body parameters that contribute to the Fourier 

transform of </> qi q 2 q 3 at large relative distances, where m' = 
m (if m is even), ml = m — 1 (if m is odd), and "round" 
rounds to the nearest integer. 



C. Asymptotics of 0^ (1-1^3) at large relative distances 



^ ( ri r 2 r 3 ) ee J ^(^a &k, exp ( g ik; • r;) , 

where Yli=i = 0- We consider the asymptotic expansions 
of <// 3 )(rir 2 r3) in two different limits: 1) the distance r be- 
tween two bosons is fixed, but the distance between their cen- 
ter of mass and the third boson, R, is large, or 2) all three 
interparticle distances, si, s 2 , and S3, are large but their ratio 
is fixed. These two expansions are respectively 



^ 3 )(r/2,-r/2,R) 



2a 2wa 2 4wa 3 24\/3W(t - 3/2) - £1 32V3wa 5 (6r - 11) - 3tt£ 2 
~ ~R + ~kR 2 ~ ~ 'ttR 3 ' + ^R4 + 2 TT 2 R 5 



C 3 (12t-25) + 12& , C4(60t- 137) + 30^4 



0(r) 



+ 



48 v / 3wa 5 (60r - 137) - 45< 2 



2ir 2 R 7 



irR 7 

27wa 2 270wa 3 



3wa 2 



18wa 3 
ttR 5 



ri? 6 



tR 7 



15a 40(2tt - 3\/3)a 2 15wa 3 



12V3w;a 5 (210r- 457) 105£ 4 



(d) 



7^ 2 i? 7 



27Ti? 7 



R 



r) + 



TTi? 4 

3(r)+ 2R 3 ' vri? 4 
180(2tt - 3y/3)a 2 315wa 3 



18V3wa 4 (12T - 25) - 9£i 
7^ 



+ 



24^; 



(d) 



945a 144(76tt - 135V3)a 2 
8#5~ + tTrs 



TTi? 6 

945wa 3 ~ 



47Ti? 7 



2irR 7 
135135a 



V R v ; 32i? 7 



TTi? 5 

4 d) (r) 



7T 2 i? 6 



(r)+0(iT M ), (45a) 



0( 3 )(n,r 2 ,r 3 ) = 1 



r 3 



^ a + 4fl2 ^ 2wft3 + 8V3wa 4 (t - 1 - 0; cot 20;) 



Si TtRiSi irB 2 Si 



r 2 B 4 



V3D 
8tt 3 B 4 



M 

i=l k 



f36wa 5 [(2i- 3) sin 30; - 20; cos 30*] - \/3a£isni30, 96wa 6 [30, 2 sin 40, + (6i - 11)0, cos 40/ 



+ 



7r 2 S 5 sin20 i 

3\/3a£ 2 0,cos40, 45V3aa d (240 l - 8 sin 40, + sin 80;) 
ir 2 B 6 sin 20, + 8ttB s sin 3 20; 

12\Z3a£ 3 sin50; " /; ' J 



P 2 (R,-s J ) + 



7r 3 B 6 sin 20 4 
V^aCs [120, cos 50, + (25 - 12t) sin 50/ 



7r 2 B 7 sin 20, 



7r 2 J B 7 sin20 



56/, / 9 N a 2 a d sin 2 0,(9 + 10 cos 20; + 2 cos 40;) ~ 1 , Rs 

f - 45 ( 2V ^ - ^) Sw^ ^ P2 ^ ' § 4 + ° {B 1 (45b) 
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where r = Ln(=yJ2/|a|), t = \n(^B/\a\), 7 = = 
1.78107 • • • , 7 is Euler's constant, and 

si=r 2 -r 3 , s 2 = r 3 -ri, s 3 = ri-r 2 , (46) 
Ri = ri — (r 2 + r 3 )/2, and similarly for R 2 and R 3 , 



B = y(sj + 4 + s|)/2, 0i = arctan(2 J R 4 /V3s J ). 
C 4 = 42 V3 wa e r s /5n - (48/5tt 2 + 32/bV3 n) wa 7 , (47a) 

£4 = [(1/10tt 3 + l/15V3 7r 2 )a - 7\/3 r s /807r 2 ] a 2 D 

- (536/25tt 2 + 124/27- 86/25^3 7r)wa 7 
+ (65/3 + 309\/3/25tt) wa 6 r s - 6wa 5 r 2 /5 

- (60tt - 87V3)a 2 a d - 9wa 4 r;/20, (47b) 



V3aL»/287r 2 - (416/21 - 8384^3/24571") too 5 
- 3wa 4 r s /7 - 3wa 3 a d r d /2. (47c) 



Ri = B sin ^j, Si = ^j-B cos and cos = 0. 

Equation (45b) derives from (40b). 

The terms up to the order R~ 6 on the right hand side 
of Eq. (45a) derive from (40a); those of the order Rr 1 , 
however, are inferred from Eq. (45b) and the continuity of 
</>( 3 )(r/2, — r/2,R) across two connected regions: r <~ 0(r e ) 
and r <~ O(R) (they join at r e « r « R), and also the 
Schrodinger equation at large R 



(ff - 3Vy4)0< 3 > (r/2, -r/2, R) = 0, 



(48) 



where H is the coordinate representation of the hamiltonian 
H for two-body relative motion [Eq. (7)]. 



III. LOW-ENERGY EFFECTIVE INTERACTION OF 3 
HARD-SPHERE BOSONS 

In this section, we consider the hard-sphere (HS) interac- 
tion. The 2-body potential V(r) = (r > 1), V(r) = +00 
(r < 1) (a = 1 in this section), and there is no 3-body po- 
tential. We numerically solve Schrodinger equation at zero 
energy in the coordinate representation, in conjunction with 
Eq. (45a) at large R, to determine D.[19] 

For 3-body configurations excluded by the repulsive inter- 
action, <j)^ vanishes; for allowed configurations, (f>^ satisfies 
the free Schrodinger equation. So J2^=i ^1 ( r i> r 2> r 3) i s 
nonzero on the boundary B between these two regions only 
(in fact, it has a <5-function singularity on B). Thus 

<t>^ (r ir2 r 3 ) = 1 - [ dc' r F(R', d)R' 2 dR' 

K J-l JRminiC) 

x {g[si,s 2 , s 3 ; 1, s-(R', c'),s+(R', c')] 
+g[s 1 ,S2,s 3 ;s-{R',c'),l,s+(R',c')] 
+G[si,s 2 , s 3 ; s_(R', c'),s+(R', c'), l] }, (49) 



for some function F(R',c'), where s,'s are defined in 
Eq. (46), d is the cosine of the angle formed by the line con- 
necting two bosons with distance 1 and the line connecting 
their center-of-mass to the third boson, 



i?mi„(c') = (|c'| + v / ^ 2 + 3)/2, 
s±(R',d) = y/l/4±R'd + R' 2 , 



(50a) 
(50b) 



and 



g(s lS2 s 3 ; s' lS2 s 3 ) = { [(s 2 + s 2 + s 2 + sf + s' 2 + s' 2 ) 2 

2 ^5 S S -y ^1^2 ^1^3 ^^2^2 ^2^1 ^2^3 ^^3^3 

— s 2 s' 2 — s 2 s' 2 2 ) — 36 (2s 2 s 2 + 2s 2 s| + 2s 2 s 2 — sf — s 2 — s 3 ) 
x (2s 1 s 2 + 2s 2 s 3 + 2s 3 s ± — Sj — s 2 — s 3 ) | 

(51a) 



is the translationally and rotationally invariant Green function, 
satisfying [Sj's are defined in Eq. (46)] 



V 2 £(si S2S3 ; s[s' 2 s' 3 ) = - J II 
i=i ^V^s 1 s 2 s 3i=1 

(51b) 

Because </>( 3 ) (rir 2 r 3 ) vanishes on £>, the unknown function 
c') in Eq. (49) satisfies an integral equation: 

(3) (rir 2 r 3 ) = at *i = l,s 2 , 3 - s T (R,c), (52) 

where — 1 < c < 1 and R > i? m in(c). 
From Eqs. (49) and (51b), we get 

(V 2 +3Vy4)0( 3 )(r/2, -r/2, R) = c)5(r- 1) (53) 

at |R ± r/2| > 1, where c = R f . 

From Eq. (53), it is clear that F(R, -c) = F(R, c). 

The 2-boson special functions satisfy Eqs. (8) at r > 1 and 
vanish at r < 1, so 

o=l, r s = 2/3, r' s = 8/15, a d = 1/45, r d = -150/7, 
a ff = 1/99225, a; = 1/1404728325. 

Applying (V 2 + 3V|/4) to Eq. (45a), and comparing the re- 
sult with Eq. (53), we get the asymptotic expansion of F(R, c) 
at large R: 

F(R,c) = FW(R,c) + QF<?\R,c) + 0(R- 8 ), (54a) 



8 



4w /24V3™ „ 12 
+ \ \nR+ — 



14 



+ 



R irR 2 irR 3 

4(72 + 43V3tt)w; , „ 10885 48^3 105 
lni?H — + 



V3tt 



/96V3w , „ 460 48 161 \ „ 
V 7T 2 9 7T 2 V3tt/ 



5005 



7T° 

27733 



54 ' 7T 3 7T 2 8a/3tt 



64(7V3tt-9)w 



7T" 



lni? + 



9 ' 7T 2 

61279 96V3 492 

135 7T 3 7T 2 



R- 7 + 



20^71" 

15865 4115 2220 



+ 



4V3 7T 



5 /80 40 V3 

2R3 + it ~ / 



)r- 



5w / 730 820 360\ 6 



120V3ii 



InR 



IT* 



R- 7 \P2(C) 



9 ^3648 1296V3\ 1 



8R 5 V 35 



9w 



7tt / i? 6 47Ti? 7 



, x 13P 6 (c) 
P 4 (c)-^, (54b) 



Pi 2) (i?, C ) = - 



6^ 24V3 /72 , 43 V3 
^R 1 ~ ?ri? 5 



-(- 



2 + 

7T / 



)r- 



+ 



144 112V3 30^^2(0) 



+ 



R-\ 



(54c) 



where Q is directly related to D: 

D = 48tt 2 Q - 192ttu;(1 - 7) - 12tt 2 

= 473.741 Oil 252 289Q - 744.946 799 290 500. (55) 

It can be shown from Eqs. (52) (at large R) and (54) that 

Q= lim / dc' / F(R', c)R l2 dR' 

- [i? 3 /3 - R 2 + 2wR/n - (4w/tt) In R] . (56) 

Solving Eqs. (52) and (56) for F(R, c) and Q numerically, 
[using {R 1 , d) + QF a (2) (R 1 , d) to approximate F(R', d) 
for R' greater than some sufficiently large value in these two 
equations, and discretizing F(R',c') for R' less than this 
value], the author finds 

Q = 5.290 844 ±0.000 005. 

Substituting this result into (55), we obtain Eq. (4). 



IV. GROUND STATE OF 3 BOSONS IN A PERIODIC 
CUBIC VOLUME 

In this section we return to the general interactions consid- 
ered in Sec. II. The 3 bosons are now placed in a large periodic 
cubic volume of side L, and 



l/L 



(57) 



is the small parameter. 



A. Ground state wave function and energy 



Let A ni „ 2n3 (X^i=i n « = 0) be proportional to the prob- 
ability amplitude that the bosons have momenta 27reiij (i — 
1, 2, 3). We shall call an integral vector, since its Cartesian 
components (along the sides of the cubic volume) are integers. 



Schrodinger equation (ft = mboson = 1) 
[27r 2 e 2 (n 2 + n 2 + n 2 ) - E] A ni „ 2n3 

g£6 ^27r€ni ) 27ren2,27ren3,27r€iij,27r€ii2 ) 27reiig-^-njn2n: 

1 

+ 2L 



3 £^27reni,27ren 2 ,27ren',27ren"^4n'ii"n3 + ■ ■ ■ — 0, (58) 



where summation over dumb momenta is implicit, and "• • • " 
stands for two other similar pairwise interaction terms, is 
rewritten as 

+ ^/ e (ki)7 e (k 2 )^^ ^kikskakik.k^kik^k;, 



+ 



^(ko^t/k, 



k 2 k'k"^'k"k 3 H = 0, (59) 



V'kjksk., = A ^2n 3 J|(27re) 3 (5(k 4 - 27ren 4 ), (60) 

nin 2 i=l 

where ELi * - 0, / k< = / , J k , k , = f f gfrgg. 



and 



J £ (k) = ^(27re) 3 (5(k - 27ren). 

n 

J e (k) = ^(27re) 3 (5(k - 27ren). 



(2tt) 3 (2tt) 3 ' 

(61) 
(62) 



In the following, k's will be taken as independent of e but 
q's be small momenta such that q/e's are independent of e. 
Also, le(k) and J e (k) will be simply written as J(k) and 
J(k), respectively. 

Let 



SO V> qi c 



A),o,o = L e , 
nL 1 (2^) 3 ( 5(q l ) + o(^ 6 )- 



(63) 
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The following expansions will be found: 

V'q,-q/2+k,-q/2-k 



k 2 k 3 1 



(s)q 
k ' 



E4 

s>-3 

? /'qiq2q3 = E ^qicuqa' 

s>-6 



(64a) 
(64b) 
(64c) 
(64d) 



s>3 



where Tl^^, S k s)q , T q ( ^ 2q3 , and £( s ) scale with e like e s 
(no/ excluding e s In e). Equation (64a) is understood in the 
following sense: we expand the Fourier transform of V'kikg.kji 
within a large but fixed spatial region in powers of e, and then 
transform the result back to the k^-space term by term, to ob- 
tain J2 S ^kik 2 k ■ Equation (64b) is similar: we expand the 
partial Fourier transform, J u V'q,-q/2+k,-q/2-k exp(ik • r), 
within a large but fixed region of the r-space, for fixed q/e, in 
powers of e, and then transform the result back to the k-space 
term by term, to obtain ^ s <S k s ^ q . 
Obviously 



= (2^) 3 <5(q 1 )(2^) 3 ( 5(q 2 ). 



(65) 



Therefore % 



T 



-6) 
kik 2 k 3 



q,-q/2+k,-q/2-k 



(2 7 r) 3 ,5(k 1 )(27r) 3 ( 5(k 2 ) and 

(27r) 3 <5(q)(27r) 3 (5(k), indicating 

Is) (s)a 
that the minimum values of s for 7?- klk2k3 and <S k are and 

—3, respectively. Since the energy is dominated by pairwise 

mean-field interactions, E <~ e 3 . 

The Fourier transform of I(k) becomes 5(r) within any 

fixed spatial region, when e is sufficiently small. Using the 

prescription in Appendix B, we thus get: 

J(k) = l + 0(e Sc ) for any s c . 

Using this result, and also noting that the interactions are 
finite-ranged, we get three specialized forms of Eq. (59) (ac- 
curate to any finite order in e): 

1 If 

-(fc 2 + fc 2 + fc 3 )7?. klk2k 3 + - / ^k 1 k 2 k 3 k' 1 k 2 k^k' 1 k 2 k?, 

1 Jk' k' 



2 7k' 



C^kik 2 k'k"^k3k'k" 



£"^k!k 2 k 3 



(66a) 



f 7 -q/2+k,-q/2-k,q,k' 1 k 2 k 3 '7<!-k 



+ 



6 j k ' k ' 
" k i k 2 

O / ^-q/2+k,q,k'k"^-q/2-k,k'k" + (q 
Z Jk> 



-q) 



(67) 



where K, S, and T are the sums of K<- s \ S<- s \ and T< s ), 
respectively, and in Eq. (66c), 1' = 2, 2' = 3, 3' = 1, and pi 
is defined in Eq. (19). 

The following statement is now true at s\ = —6: 
Statement si: The functions 7q\ s q2q3 (for s < si), 5 k s ^ q 

(for s < si + 2), ^ k s i ) kak3 (for s < Sl + 5), and S( s ) (for 
s < si + 8), have all been formally determined. 
Now do the following expansions for —6 < s < s\: 



q,-q/2+k,-q/2-k 
-'kikaka 



t(n,s—n) 
k i 



E./l(n,s—n) 



ik 2 k 3 



(68a) 
(68b) 



where q/e's and k's are indepedent of e, ki <~ fc, and 

_ £ n fcs -n f"("' s "") ^ e ™fc s -«. For the same re a- 

q,k ' Kik 2 K 3 

son as Eq. (25a), we obtain the following asymptotic expan- 
sions at small k's: 



j(«l+3)q 
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(si+6) 
kik 2 k 3 



-3 

X 'q.k 

=-si-9 

-6 

y / ' '"- 1 k .; k 

=-si-12 



^ +3 ' m) +0(e 5l+3 fc- 2 ), (69a) 
^t 1 „t 6 ' m) +0(e 5l+6 fc- 5 ). (69b) 



Extracting all the terms that scale like e s i+ 6 from Eq. (66a), 
we get 



fcf 



K 3J /C k!k 2 k 3 



kikaksk'jkfjks /v k', k^k^ 



(ai+6) 



+ 2 / ^kik 2 k'k''^ k3 



(si+6) 



k 3 k'k" ~r 



^ /C k!k 2 k 3 ' 

(70) 



3<n<si+6 



where the right hand side is already known. This equation and 
Eq. (69b) [which expands T^^ 3 to the order fc~ 6 at small 

k's] are sufficient to determine T^^ 2 ^ 3 , due to the 3-body 
version of the uniqueness theorem [analogous to Eq. (13)]. 

The above information for 6> k and 7£. klk2k3 is sufficient to 
determine each term in Eq. (66b), except the first term, to the 
order e Sl+3 [note that 7(q) ~ e°], so (i7S (si+3)q ) k is known; 
taking into account Eq. (69a), one can determine 6> k 1+3 ^ q . 
The sum of all the interaction terms in Eq. (66c) can now 



(HS q ) k + (3g 2 / 4 - E)S£ + I(q)W q = 0, (66b) be determined to the order e Sl + 3 , with a result 



-qiq 2 q 3 



+ 



\{(il + <£ + ql)-E 



7q iq2q 3 

A »=i Jv 

+ ^/(qi)/(q2) / ^ f7q 1 q 2 q3k' 1 k 2 k3^-k' 1 k 2 k 3 = 0, (66c) 
Jkjk 2 



c o]li=i( 27re ) S(<h), where c' qiq2q3 vanishes at qi = q 2 = 
0. Both c qiq2q3 and Co are known up to the order e Sl+3 . 

Because of Eq. (63), = T qiq2q3 -nLi(2^) 3 <5(qO 

vanishes at qi = q 2 = 0. Solving Eq. (66c), we get E = c e 6 
(so now E is known to the order e Sl+9 ) and 



qiq2q 3 



G a 



qiq2q 3 



"qiq2q 3 



). 



(71) 



10 



where G qiq2q:3 = 2/(qf + q\ + g§) ~ e~ 2 . Because E ~ e 3 , 
^ iq2q3 ~ e 2 , a; > -6, and £ and ^ iq2q3 are known to the 
orders e Sl+9 and e Sl , respectively, we know £^ iq2 q 3 to the 
order e Sl+3 . So from Eq. (71), we can now determine ^ iq2 q 3 
(and thus ^q iq2 q 3 ) to the order e Sl+1 . 

The truth of Statement (s\ + 1) is now established. 

Repeating the above routine, one can formally determine 
^4nin 2 n 3 and E to any orders in e. The following are the step- 
by-step results of this program. 

Step la. TZ (0 \ TZ (1 \ and TZ^ satisfy the zero-energy 
Schrodinger equation [because of Eq. (70)]. They will all be 
determined in this paper. From Eq. (65) we get i^k = 
II; ,;2rr ! V;k,,.So 

< ) k 2 k3=^ 3 1 Us- (72) 

Step lb. t'^'-V = (27r) 3 (5(q)(2 7 r) 3 ( 5(k); (HS^) k = 
0. So 



Step 3a. From t'-^^ and Eq. (70), we get 



<Us = -( 3 «w*)< kak8 - 



(3) 



(79) 



Step 3b. From (HS^ 2 ^ = 0, and t' ( ~ 2, 3) , we get 



4" 3)q - (27r) 3 <5(q)</»k. 



(73) 



Step 2. With the help of Eqs. (12), and noting that 

Ei=i = °> we s et 



E = 12nae 3 + o(e 3 ) 



(74) 



4 = [-(aiaeA)(2^) 3 ( 5(q)-87r a J(q)/ (Z 2 ]0k. (80) 
Step 4. This is similar to Step 2. 

E = 127rae 3 (l-aiae/7r) + 0(e 5 ), (81) 



T(-4) —7 r 1 \ ^ 9 2 / 2 4 2 

i=l 

3 

+ (a 2 /7r 2 e) ^(a im r 2 + m 7 4 ) J(p 4 )(2^) 3 <5(q 4 ), (82) 



n' li27 t0,n' 1 +n^0*=l 

The three subscripts of J qi q 2 q 3 are always subject to the 
constraint ^ i=1 qi = 0. In Eq. (82) and in the following, 

G„ in2 „ 3 = 2/{n\ +nl+ nf), and 



T ^L=Y,-^ na / p ^ J ^ 2n ^ 5 ^- (75) n, = q,/2vre, n = q/2vr e , m, = Pi /27re. (84) 



Using the method of Appendix B, we get 



T 



( q/2 + k,-q/2-k,q = [ " 8™/^ " ™(q ' V fe ) 2 



+ 0(e 2 )] (2^) 3 <5(k) + (27r) 3 <5(q) [ - 47ra/fc 2 

- ai ae(27r) 3 (5(k)/7r + 27rae 3 V^(27r) 3 (5(k)/3 + 0(e 5 ) 



-'q,-q/2+k,-q/2-k - Z^=-2 f q,k ' Wnere 

4k 1 '" 35 = (« 2 A^){(«? + a2 )(2 7 r e ) 3 ( 5(q)(2 7 r) 3 ( 5(k) 

+ J(q) [2p A1 (n) + 2 ai /n 2 -6/n 4 ] (27r) 3 <5(k)} , (85a) 
*q ( °k~ 4) = [l6^V(2^) 3 <5(q)+40^ 2 a 2 J(q)]Z/fc 4 , (85b) 



(76a) 



(85c) 



(-5) 



kik 2 k 3 



X-(4W«» 2 )(2 7 r) 3 5(k i )] 



i=i 



- (3aiae/7r)(2 7 r) 3 (5(k 1 )(27r) 3 ( 5(k 2 ) + 0(e 3 ), (76b) 



The actual formulas for (s = —4,-3,-2; see Steps 6 
and 8 below) are not needed in this paper. For s > 1, 

p As (n) = 20 Aa (n) + ^(n)/n 2s , (86a) 
9 As (n) = [(™ 2 + m • n + n 2 )m 2s ] ~\ (86b) 

m^0 



h = (k 2 -k 3 )/2, and similarly for 1 2 ,1 3 . (77) ^(n) = lim V (m 2 + m • n + n 2 )" 1 - 4tt7V. 



|m|<JV 



}3„ 



a s = Km V — - I d*n— (s < 3/2), (78a) 



(86c) 



tl . 5 =lim[( ]T |n|- 3 )-47rlniV], (78b) 

TV— >oo ^ — ' 

(78c) 



^ ( rk 4 2 ) k3=^ lfe2fe3 E 327rV A t 2 



n/0:n<JV 

a s ^J2\^ 2S (*>3/2). 



i=i 



n^0 

The above lattice sums are evaluated in Appendix C. One can 
easily extract iq^ ° ^ and i k lk 2 ~k 3 ~'^ fr° m Eqs. (76). 



+ 12 ai a 2 e]T(2^) 3 <5(k 4 )// 2 



+ 3ir- 2 (2(3 1A + a 2 - 3a 2 )a 2 e 2 (2 7 r) 6 ( 5(k 1 ) ( 5(k 2 ) + 0(e 3 ), 

(87) 
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from which one can easily extract ik|k 2 k 3 ^ ■ ^ ere 

[3 1A = lim V c- r i n W A {n)/n 2 + aVZtt 3 /t 1 2 (88) 

is evaluated in Eq. (C18). 

Step 5a. From t'^~^ and Eq. (70), we get 

<U 3 = 37r- 2 (2/3^ + a? - Sa^aV^ 3 ^. (89) 

Step 5b. From (if^" 1 )^ = 0, and t'^ 1 '®, we get 

4" 1)q = (a 2 /TT 2 e){(al+a 2 )(2TTe) 3 5(<D 

+ [2p A1 (n) + 2 ai /n 2 - 6/n 4 ] J(q)}0 k . (90) 

Step 6. This is similar to Steps 2 and 4. 

E = l2-Kae 3 [l-a l ae/-K+{a\+a 2 )a 2 e 2 /-K 2 ]+0{e & ), (91) 

3 

^qiq2q3 = Z 71 " 6 )^qiq2q3Cnin 2 Ti3 n i 

»=1 

- (2a 3 /7T 3 e 3 ) Jq 1 q 2 q 3 G„ 1 „ 2ll3 [p A l( n ») + ~ 3 /^] 

+ { [-Ap Al {m l )/m 2 -{a 2 +a 2 )/m 2 + 2 ai /m 4 + 15/mf] 
x a 3 /7r 3 + «o} J(p 4 )(2^) 3 <5(q 4 ). (92) 

In preparation for the subsequent steps, we derive 

t >V,-S) = { _ 327r 2 wa3Zi/iai{k)/k 3 

+ [2u + 16wa 3 ln(27re|a|)](27r) 3 (5(k) + (2a 3 /tt 3 ) 

x [ - PAAi( n ) - a iPAi( n ) + 3p A2 (n) + 3p B i( n ) 
- (a 2 + a 2 )/n 2 + 6 ai /n 4 - 9/n 6 ] (2vr) 3 (5(k)} J(q) 
+ { - 327r 2 wa 3 Z 1/ | Q |(fc)/fc 3 + [I6wa 3 ln(27re|a|) - w 
+ (15a 3 + a x a 2 - a\ - 4a 1A1 )a 3 /tt 3 ] (27r) 3 <5(k)} 
x (27re) 3 5(q), (93a) 

ik" 4 - 47k ^(q) - (967r aia 3 e Z/fc 4 )(2vre) 3 5(q), (93b) 

where p B1 (n), p AA1 (n), and a 1A1 are defined as follows. 
For s > 0, 

Ps » = 20 Bs (n) + ^ B (n)n- 2s , (94a) 
Bs (n) = [(«* 2 + m • n + n 2 ) 2 m 2s ] "\ (94b) 

m^O 

VF B (n) ee ^(m 2 +m-n + n 2 )- 2 . (94c) 

all m 

PAAi( n )= J im [ 2 y] (m 2 + m- n + n 2 )"V^i(m) 
- 87r 3 wln7V] + W A (n)p A i(n). (95) 



It can be shown that at large n 

p A1 (n) = 7T 2 w/n + 2 ai /n 2 + 0(n~ 4 ), 
p AA1 (n) = -87r 3 w;lnn - ir 3 w(7ir/V3 - 8) + 2ir 2 wai/n 
+ 0(n- 2 ). 

We define some 2-fold lattice sums: 

a 1A1 =^^{1 J2 n- 2 p A1 (n)}-47r 3 w\nN}, (96a) 

n^0;n<A r 

a sAs , EE^n- 2 V^(n) (s + s'>3), (96b) 

n^O 

« sBs , EE^n- 2s p B v(n) (s + s'>2). (96c) 

n^O 

a 1A1 , a iA2 , a 2A i, and aisi are evaluated in Appendix C. 
Step 7. From 

(ff5<°^) k + (-6W k + ^ 0) ) j(q) 

+ (-127ra0 k + TV^ 0) )(27re) 3 (5(q) = (all k), 

4 0)q = £ Ck s) + °( e °^ 2 ) ( smallk )> 

s=-5 

where is defined in Eq. (22b), we get 

5 k ° )q = {[l6wa 3 \n(e\a\) + 2a 3 p AA1 (n)/ir 3 ]cf> ]i + d k 
+ IOtto^} J(q) + { [I6wa 3 ln(e|a|) + C a 3 + 3na 2 r s 
- 3u Q ]fa + d k + 6™./k}(2^) 3 <5(q), (97) 

where d k is defined in Eqs. (36), 

C = 16wln(27r) + (15a 3 + a x a 2 -a\- Aa 1A1 )/Ti 3 
- (14tt/\/3 - 16) 
= 95.852 723 604 821 230 29, (98) 

and 

p AA1 (n) = - p AA1 (n) - tt 3 w (7ir/V3 - 8) + 8tt 3 u; ln(27r) 

- a iPAiH + 3p A2 ( n ) + 3p B i(n) 

- (a 2 + a 2 )/n 2 + 6ai/n 4 -9/n 6 . (99) 

See Appendix C for the evaluation of Co (and C\ in Step 9.). 

Step 8. This is similar to Steps 2, 4, and 6. With the help of 
Eqs. (12) and (39), we get 

E = 127rae 3 {l - aiae/n + (a 2 + a 2 )a 2 e 2 /7r 2 

+[l6wa 3 ln(e|a|) + C a 3 + 3Tra 2 r s ]e 3 }+De 6 + 0(e 7 ), 

(100) 
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T £L = E [J(p,)(2 7 r) 3 ( 5(q,)T(-^(m 4 ) 
i=l 

+ -Wqs EL G n 1 n 2 n 3 ^ ( " 2) (n,)] , (101a) 

T^ _25 ^(nij) = — aiauoe/7r — auoe/irm 2 — De/47r 2 TO 2 

- 3a 3 r s e/m 2 + { - [Ap AAl {m. t ) + A8n 3 w ln(e|o|) 

+ 7r 3 C ] /m? - [I2p ill (m i ) + 9a? + 6a 2 ] /m 4 

+ 3ai/TO^ + 45/TOf}a 4 e/7r 4 , (101b) 

T ( " 2) (n 4 ) = -2au /ne 2 nl (101c) 

^"^(m) = -2[p AA1 (m) +^w\n{e\a\)]a A /^e 2 

-D/12n 2 e 2 , (lOld) 

T 2 ( - 2) (n 4 ) = -6[p A1 (n t ) + - 3/n 4 ]a 4 A 4 e 2 , 

(lOle) 

T 3 ( - 2) ( ni ) = 18a 4 /7r 4 e 2 n 2 . (lOlf) 
In preparation for the next step, we derive 

<ik~ 3) = ^7k^(q) + {967r U ; ai a 4 Z 1/ | a |(fc)/fc 3 
+ [ - (Wwajir)^ ln(e|a|) + C[a 4 - aiL>/47r 2 

- 3a ia 3 r s ] (27r) 3 ,5(k)}e(27re) 3 (5(q), (102) 
C( = Ti-- 4 [ - 45 1AA1 - 487r 3 wa 1 ln(27r) - Tr 3 ai C 

— 12a 2A1 — 9a 2 a 2 — 6a 2 + 3a 1 a 3 + 45a 4 ] , (103) 

where 

= { [ E m_2 ?W m )] 

m^0:m<AT 

- 327rW[ln(27riV) - l] + 12n 3 wa 1 lnN} (104) 

is a 3-fold lattice sum (evaluated in Appendix C). 
Step 9. From 

(FS (1)q ) k + [24a ia 2 e 4 k - 3a 1 ae 4 W^ >) /tt] (2^) 3 <5(q) 
+ { 6 [n 2 p A1 (n) + + n" 2 ] a 2 e^ k - 3aiae W^ 0) /tt } J(q) = 

and [note that ^" 7) = ^" 6) = 0] 

4 1)q = E & + O^k- 2 ) (for small k), 

s=-5 

we get 

5 k 1)q = { [ - ^feue^/TOa 4 ln(e|a|) + da 4 - ai D/4TT 2 
— 6ctia 3 r s + 6otiau a /ir] </> k — 6aia 2 / k — 3aiad k /n} 
xe(2^) 3 ( 5(q)+y(q,k)J(q), (105) 

where 

Ci = C7J +6tuai(7/v^-8/7r) 

= 810.053 286 803 649 420, (106) 



and the actual formula for Y (q, k) is not needed in this paper. 

Step 10. Substituting the latest results for <S q and 7\L qiq2q3 
into Eq. (66c), extracting all the terms that contain the factor 
(27r) 6 <5(qi)<5(q2) from this equation, noting that 7q iq2q3 — 
(27r) 6 <5(qi)(5(q 2 ) vanishes at qi = q 2 = 0, and using 
Eqs. (12) and (39) to simplify the result, we get 

E = 127rae 3 {l - a^e/ir + {a\ + a 2 )a 2 e 2 /iT 2 
+ [l6wa 3 ln(e|a|) + C a 3 + 3ira 2 r s ] e 3 
+ [ - (96wai/ir)a 4 ln(e|a|) + da 4 - 6aia 3 r s ] e 4 } 
+ (e 6 -6a 1 ae 7 / 7 r)L) + 0(e 8 ). (107) 

Substituting the numerical values of the lattice sums (see Ap- 
pendix C) to the above equation, we get Eq. (5). 

We will not proceed to determine Tq~q]q 3 in this paper. 



B. Results for Wu's parameter £3 

Wu computed the three-boson energy in the large periodic 
cubic volume to order L~ 6 , and he left an unknown parameter 
£ 3 for the three-boson interaction strength at low energy [4]. 
£3 is needed to determine the full order-p 2 correction to the 
many-body energy [4] . 

Comparing our Eq. (100) with Wu's result, Eq. (5.29) of 
Ref. [4] (note that the unit of mass in [4] is 2m boS on = 1), 
we find that Huang and Yang's constant [20] C = —ai/ir — 
2.837 • • • , [21] different from the number 2.37 which was first 
provided in Ref. [20] and then adopted by [4] ; the symbols £ s 
in Ref. [4] equal a s /ir s ; the symbol £3 in Ref. [4] is now 
expressed in terms of D of the present paper: 

£3 = D/Utto 4 + 3irr s /a 

- 4a 1A1 /ir 3 + 16wln(27r) - (Uw/Vs - 16)w 

= D/l2Tra 4 + 3irr s /a + 73.699 808 371 935 4035. 

(108) 

For hard-sphere bosons, D is given by Eq. (4) and r s = 

2a/ 3, so 

£3 = 126.709 37 ± 0.000 06 (for hard spheres). (109) 



C. Results for the ground state wave function 

There are 6 expansion formulas for A ni „ 2n3 , each of which 
is valid in a sub-region of the discrete momentum configura- 
tion space: 

A oao = L 6 , (110a) 
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A ,n,-n = -aL 5 /im 2 + (ai/n 2 + 1 / n 4 ) a 2 L 4 / 'tt 2 + { - [a\ + a 2 + 4p A1 (n)] /n 2 + 2 ai /n 4 + lb/n 6 }a 3 L 3 /n 3 + u L 3 
+ {[- 4 pW n )/™ 2 + ^ 3 wn- 2 ln(L/\a\) - 7r 3 C /n 2 - (9a 2 + 6a 2 )/n 4 - l2p A1 (n)/n 4 + 3ai/n 6 + 48/n 8 ]a 4 /7r 4 

- (ai + n- 2 )au /7r - D/ATt 2 n 2 - 3a 3 r s /n 2 }L 2 + OIL 1 ), (110b) 

3 

A nin2n3 = ]T{2a 2 L 4 /7r 2 n 2 -2[^>0+^ 2 -3/n 4 ]^ 

x G ninans + {6a 3 i 3 /7r 3 n 2 - 6[p A1 (n 4 ) + 2 ai /n 2 - 3/n 4 ] a 4 i 2 /7r 4 }G 2 in2 „ 3 + l8a 4 L 2 G 3 nin2n Jir 4 n 2 - 2au a L 2 /i,n 2 

+ 0(L 1 ), (110c) 

A),n,-n = {L 3 -a 1 aL 2 /n+(a 2 1 +a 2 )a 2 L/TT 2 -16wa 3 \n(L/\a\)+C a 3 +3na 2 r s -3u +[(m^ 

- ai D/4ir 2 - 6a ia 3 r s + Gaxauo/ir]!,- 1 }^ + (l - 3aia/7rL)d k + (67ra - 6a ia 2 /L)f k + 0(L~ 2 ), (HOd) 



^n,-n/2+N,-n/2-N = { - 2aL 2 / W + [2p A1 (n) + 2ai/n 2 - 6/n 4 ]a 2 L/i: 2 - 16«;a 3 \n(L/\a\) + 2a 3 i o 4A1 (n)/^ 3 }</> k 

+ 107ra4k +dk + 0(L- 1 ) 7 (HOe) 

A Ni n 2 n 3 = [1 - 3a ia /nL + 3(20 1A + a\ - 3a 2 )a 2 / 'n 2 X 2 ] + O^ 3 ), (1 lOf) 



where n's are nonzero vectors of order unity, N's are large 
vectors of order L / max(r e , | a | ) (r e is the range of the interac- 
tion), k = 2ttN/L, ki ee 2nNi/L, and G ni „ 2 „ 3 ee 2/{n\ + 
n 2 + n D- T ne formulas for A 0j o,o. A) ; 11 . 11 ' 

and A nin2ll3 

are extracted from the above results for ^q iq2q3 ; those for 
A),n,-n and A n ._ n / 2 +N.-n/2-N are extracted from S^; the 
formula for ^4niN 2 n 3 results from the expansion of 7^kik 2 k 3 - 
The numerical constants in the above formulas are computed 
in Appendix C. 

D. Momentum distribution 

The expectation value of the number of bosons with mo- 
mentum 2im/L is Af n = c^ alln , l^n^/^n-n/l 2 for some 
constant c such that ^ all n 7V n = 3. 

For any nonzero integral vector n of order unity and any 
large integral vector N of order L/ max(r e , |a|), we have 

Wo = c\A 0a0 \ 2 +C ^ l^0,n',-n'| 2 +C ^ I-Vn',-N'| 2 , 

0<n'<JV c N'>N C 

(111a) 

7V n = 2c|A ,„ ; _ n | 2 + C ^2 IAi,-n/2+n\-n/2-n'| 2 
n'/±n/2, n'<JV c 

^n,-n/2+N',-n/2-N'l ) (Hlb) 

N>>N C 



-A/jsr = 2c| j 4 ,n,-n| 2 + 2c ^ l^n<,N,-„<-N| 2 

n'/0, n'<N c 

+ c ^2 |-4n',n -n'-n| 2 , (Ulc) 

N'>N C , |N'+N|>7V C 

where N c satisfies 1 <C N c <C Lj max(r e , |a|); inEq. (111b), 
n' ± n/2 and N' ± n/2 have integral Cartesian components. 
Substituting Eqs. (110) and using Eq. (14), we get 

Wo/c = L 12 + a 2 a 2 L 10 /ir 2 - 2{a x a 2 + a 3 )a 3 L 9 /Tt 3 
- 2nu L 9 - 2na 2 r s L 9 + 0(L 8 ), 

Wn/c = 2a 2 L w /ir 2 n 4 - A( ai /n 4 + 1 / n 6 ) a 3 L 9 / tt 3 
+ 0(L S ), 

A/n/c = 2L 6 (1 - 2a ia /nL)\^ /L \ 2 + 0(L 4 ), 
Solving J2 d \\ n -^n = 3 [using Eq. (14) again], we get 

c = 3ZT 12 [1 - 3a 2 a 2 /Tr 2 L 2 + 6( ai a 2 + a 3 )a 3 /ir 3 L 3 

+ mu /L 3 + 6ira 2 r s /L 3 + 0(iT 4 )] , (112) 

and 
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7V /3 - 1 = -2a 2 a 2 /n 2 L 2 + A{ ai a 2 + a 3 )a 3 /ir 3 L 3 + (43?u + Aira 2 r s )/L 3 + 0(iT 4 ) 

= -3.350147643 a 2 /L 2 + (43?u + A^a 2 r s - 17.926831164 a 3 )/L 3 + 0(iT 4 ), (113a) 

7V n = Qa 2 /iT 2 L 2 n A - (\2a 3 /ir 3 L 3 )( ai /n A + l/n 6 ) + 0(iT 4 ) 

= (0.607927102/n 4 )(a/L) 2 + (3.449740068/n 4 - 0.387018413/n 6 ) (a/L) 3 + 0(iT 4 ), (113b) 

A/n = 6iT 6 (l - 2a 1 a/nL)\(f >2 ^ /L \ 2 + 0{L~ & ) = 67r 6 (l + 5.674594959 a/L)\^ N/L \ 2 + 0(iT 8 ). (113c) 



One can derive higher order results for TV's from Eqs. (110). 

Equation (1 13a) shows that the depletion of the population 
of the zero momentum state depends on the parameter u at 
the next-to-leading order in the volume expansion. 

Given <f)(r), one can compute u and r s from the formulas 



u = J [</>(r) - (1 - a/r)} d 3 r, (1 14a) 

-2na 2 r s = J [\<f>(r)\ 2 - (1 - a/r) 2 ] d 3 r. (114b) 



Because <p(r) = 1— a/r outside of the range of the interaction, 
the integrands are nonzero within the range only. Because 
there exists a two-body potential V(r) = (/> _1 (r)V 2 (/>(r) for 
any function (j)(r), Eqs. (114) indicate that u is in general 
independent from a and r s . 

After repeated failure to find a general relation between 
3?uo an d an Y finite number of parameters in Eq. (9), the au- 
thor conjectures that such a relation does not exist at all [22]. 
Because the EFT [9] is formulated in terms of parameters in 
Eq. (9) along with 3-body, 4-body, ...parameters, the two- 
body parameter u is absent in the EFT. 



E. Generalization to N bosons 



tact interaction parameter 7/3 (^) in Ref. [23] satisfies 

m (\a\- 1 )-48a A (4Q + 2K)/TT 2 
= D + 12ir 2 a 3 r s - 48aV A1 /tt 2 

+ 247ru;a 4 [8 ln(27r) - (7tt/V3 - 8)] 
= D + 127rVr s + 2778.417 318 626 973 645a 4 , (115) 

where 4Q + 21Z is a number defined in Ref. [23]. It is the 
combination r/ 3 (/i) - 48a 4 (4Q + 21Z) /it 2 that appears in their 
formula [23] for E. In summary, 

E = p { i/L-a,N) jjr- L m R 



+ 



Af\ D + \2ix 2 a 3 i 

T 6 



+ 



Af\ 8TT 2 a 3 ° 



L 6 



+ 0(L- 7 ), 



(116) 



where Pisa well-determined [23] power series in 1/L with a 
and M as the only parameters, and (™) = jy^zjy, ■ 

In addition to the energy, the present author obtained the 
following results for the momentum distribution: 

x = Af /Af- 1 = -(AT -l)a 2 {a/nL) 2 
+ 2(Af- l)(nu + na 2 r s )/L 3 

+ 2{N - 1) [ ai a 2 + (27V - 5)os] {a/nL) 3 + 0(iT 4 ), 

(H7a) 



1. Results 

To understand the ground state energy and momentum dis- 
tribution of dilute Bose-Einstein condensates (BECs), the au- 
thor generalized the above calculations to J\f bosons (J\f = 
1,2,3,4...). 

Solving the TV-boson Schrodinger equation perturbatively 
in powers of 1/L, using the same finite-range interactions as 
above, the author obtained the volume expansion for E up to 
the order L~ 6 , and found that it exactly agrees with Beane, 
Detmold, and Savage's result for E [23] which was derived 
with zero-range pseudopotentials [23], if the three-boson con- 



W„ =M{M- l)n- 4 ( a /7rL) 2 

- 2N{N - 1) [ai/n 4 + (2VV - 5) /n 6 ] (a/irL) 3 
+ 0(L- 4 ), (117b) 



A/n = Af(Af - 1)L~ 6 (1 - 2a ia /irL)\0 2nN/L \ 2 + 0(iT 8 ). 

(117c) 

Equations (117), unlike E, can not be derived from the zero- 
range pseudopotentials or the effective field theory. [Even in 
Eq. (117b), the L~ 4 correction will contain the short-range 
parameter u .] 
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2. Implications for dilute Bose-Einstein condensates in the 
thermodynamic limit 

Let a > 0. Let p = Af/L 3 . At large Af there are two differ- 
ent low-density regimes (and an intermediate regime between 
the two), depending on the box size L: 

• L^> Afa, so that L is small compared to the BEC heal- 
ing length £ ~ (pa) -1 / 2 ; 

• Af 1/3 a < L < A/a, and the system is a dilute BEC 
near the thermodynamic limit (L ^> £). 

The 1/L expansions for the energy and the momentum 
distribution are valid in the first regime only. In the sec- 
ond regime, they diverge like ^2^ Ci(aAf / L) 1 . Neverthe- 
less, one can infer many properties of the BEC in the second 
regime. 

Now we tentatively take the thermodynamic limit (p fixed, 
Af and L large) of the 1/L expansions of E = E/Af and 
x [see Eq. (117a)]. Each term that remains finite is retained; 
each term that diverges must be rendered finite by a resumma- 
tion [3] that includes all similar, but higher order (and increas- 
ingly more divergent) contributions [3]. 

Thus the mean-field energy term for E is reproduced. The 
logarithmic term w — 32Trwa 4 p 2 ln(L/a) is rendered finite by 
a resummation that changes L to 0(£), yielding precisely the 
same logarithmic term as in Eq. (1). The leading nonuniversal 
terms (in the sense that parameters other than a contribute) 
become p 2 (D/6 + 2-K 2 a 3 r s ). Comparing these findings with 
Eq. (1), we infer that 

E = 2Trpa[l + (128/15V^)(pa 3 ) 1/2 + 8wpa 3 ln(pa 3 ) 



+ (L>/127ra 4 + 7rr s /a + C E )pa 3 



(118) 



plus higher order terms in density p, where C E is a universal 
constant that remains the same for all Bose gases. Because E 
was computed by Braaten and Nieto [7], a comparison will be 
made between Eq. (118) and Ref. [7] (in Sec. V D). 

The contributions to x that depend on a only are diver- 
gent in the thermodynamic limit. After a resummation that 
includes higher order terms in Afa/L (analogous to Eq. (55) 
of Ref. [3]), they must reproduce Bogoliubov's well-known 

. . . The leading nonuniversal 



formula x = 

term becomes p(2jRu + 2na 2 r s ); because this term comes 
from the short-range behavior in two-body collisions, it must 
extend into the thermodynamic limit. Combining these obser- 
vations with an EFT prediction for the condensate depletion 
[9] (which should be valid through order p 1 at u = a 2 r s = 0, 
at the very least), we infer that 



x = 



'pa 3 + p(29?u + 27ra 2 r s ) + C x pa 3 (1 19) 



plus higher order terms in p, where C x is another universal 
numerical constant. Thus the nonuniversal effect in the quan- 
tum depletion of the condensate is larger than the EFT pre- 
diction [9] by a factor of order (pa 3 ) -1 / 2 at low density, if 

1/3 

u ' ~ r s ~ a. This disagreement is entirely caused by the 



fact that the momentum distribution at k <~ l/r e S> (pa) 1 / 2 
is approximately p 2 |</>k| 2 , rather than a structureless function 
w 167r 2 a 2 p 2 /fc 4 as implied by the EFT. 

The condensate fraction (namely Afo/Af, or 1 + x) of a di- 
lute Bose gas of hard spheres (for which 2u + 2na 2 r s = 
^-a 3 ) is thus slightly greater than that of a Bose gas with 
r e <C a [for which |2u + 2ira 2 r s \ <C a 3 ] by about ^fpa 3 
at zero temperature, if the two gases have the same number 
density and scattering length. 



V. LOW ENERGY SCATTERING AMPLITUDES OF 
THREE IDENTICAL BOSONS 

In this section we compute the T-matrix elements of three 
identical bosons at low energy. We will compare our results 
with a similar calculation in Ref. [7], to establish the relation 
between the parameter D in the present paper and the three- 
body contact interaction parameter 53 (k) in [7]. 

The interactions are the same as in Sec. II. 

The T-matrix is denoted with roman type T below, since the 
italic T is already used for the wave function components. 



A. 2-boson scattering amplitude 

The two-boson T-matrix can be expressed in terms of the 
scattering phase shifts [24] (h = mboson = 1): 

T(&,b-q) = (47r/u?) ]T [e 2i4 <^ - l] (21 + l)P ; (b • q), 

(=0,2,4.... 

(120) 

where ±b and ±q are the momenta of the two bosons before 
and after the scattering, respectively, and q = b. 

At small q we apply Eq. (9) to the above formula and get 

T(6,b-q) = -87ra + i87ra 2 g + 87ra 2 (a-r s /2)a 2 

-i87ra 3 (a-r s )a 3 + 0(a 4 ). (121) 

The first two terms in this expansion agree with Ref. [7], and 
all higher order corrections disagree, since the effective range 
r s is not included in [7]. Equation (121) agrees with Ref. [9] 
where r s is taken into account. 



B. 3-boson scattering amplitude 

Equation (121) can be alternatively derived from a system- 
atic perturbative solution to the two-body Schrodinger equa- 
tion at incoming momenta ±b and energy E ~b 2 : 

+ (terms that are regular at q 2 = E), 
where — it] specifies an outgoing wave (rj — ► + ). 
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Similarly, from the stationary wave function ^q^m de- 
scribing the scattering of three bosons with incoming mo- 
menta bi, b 2 , and b 3 , and energy 

E = (bf + b\ + &§)/2, (122) 

one can extract the three-boson T-matrix elements: 

*qiq 2 q3 = E S (<& ~ bpi)<5(q 2 - b P2 ) 

P 

+ ^f l9293 [T(b 1 b 2 b 3 ;q 1 q 2 q 3 ) 

3 

+ T(/i i ,h i -p j )(2 7 r) 3 5(q, -b i ) 
+ [terms that are regular at (gf + gf + (jf) /2 = E]. (123) 

Here P refers to all 6 permutations of "123", 2~^ 3 =i b « = 
X^j = i q« = 0, Pi is defined in Eq. (19), and 

G E m = M + <& + «i)/2 - £ - _1 , (124) 

hi = (b 2 — b 3 )/2, and similarly for h 2 , h 3 . (125) 

Let bi ~ <7i ~ g be small. We determine * perturbatively 
next. 

The equations for here are formally identical with 
Eqs. (66) in Sec. IV, except for three differences: 1) the box 
size L = oo here, so I(q) = 1, 2) E ~ g 2 here, instead of 
1 / i 3 in Sec. IV, and 3) the leading contribution to ^q^^ is 

T ^L = ^jf- E S ^ ~ b P1 )<5(q 2 - b P2 ). (126) 

p 

Naturally, 5 , q 1 q 2 q 3 — * Wiicuqa m tne limit E, bi — » 0, 
and the calculation here will be a generalization of Sec. II to 
nonzero incoming momenta. So we use the same symbols T 



(^ T) and S as in Sec. II, in the asymptotic expansions 

oo 

^qiq2q3 = E/ ^qitfeqs' (127) 

s=-6 

oo 

= *q,-q/2 + k,-q/2-k = £ ^ 28 ) 

s=-3 

where Tqf q2 q 3 and 5 , ^' q scale like q s (not excluding 
g s ln m g). When E, bi — > but q's and k are fixed, the (com- 
plicated) results for T( s ) and S^' in this section will reduce 
to the much simpler ones in Sec. II. 

When the three bosons all come to a region of size <~ r e 
(radius of interaction), effects due to a nonzero E are small. 
So at momenta k t ~ l/r e > \f~E~ [25] 

*k lk2 k 3 =^ 3 1 ) k2k3 +0(Vl). (129) 
Employing the same systematic expansion method as in 
Sees. II and IV, we obtain the following results (listed in the 
same order as they were obtained): 

1 3 



i,3 = l 



5 k _2)q = - o E [ 87raG qb, +i«^(2^) 3 <5(q-b 4 )]^ k , 
J i= i 

(130c) 
(130d) 



, 3 

4 _1)q - 3 E{ [8™ 2 (i^- V3? 2 /4 - £ - ir;)Gf bi +647r 2 a 2 C f (q.b*)] &c+/i?(27r) 3 <5(q-b i ) [^-a{a-r s /2)4>^4> { ^, 



(130e) 



= ^f^E [W^a 3 (yfi$/4 - E -'vi -i^)G q 3 3bi -2567rVcf (q^b,) + 4™ 2 (a - r s /2)fe 2 (2^) 3 ,5(q, - b 4 ) 



3 



+ -f J2 ( 2 *f 6 (<h - bi), (130f) 
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4° )q = 4^E [ 6 ™ + 2 ™^ E + i^(2-) 3 ^(q - b 4 )] A + ^ E ^qUgl 

i=l i=l 

. 3 

+ 3 S [647r 2 o 3 (V3« 2 /4--B-i»? - toi)cf (q, b<) - 256^ 3 a 3 cf (q, b i; lap 1 ) + 87ra 2 (a - r s /2)(£ - 3fl 2 /4 + h 2 )G, 

i=l 

+ i 87ra 3 hi y/Zq 2 /4 - E - vqG^. - 3ira 2 r s - (Uir/VS - 16)wa 3 + ia 2 (a - r s )/i 3 (27r) 3 (5(q - b,)j </> k , (130g) 



E 

qbi 



t { - 2 ^ = n E 

qiq2q3 w 919293 



{ - D - ^ [647r 2 a 3 (^3<z 2 /4-£-i77 - i/^cf (q,-, b 4 ) - 256^ 3 a 3 cf (q,, b i; lar 1 ) 

i,j = l 

+ 87ra 2 (a - r s /2)(£ - 3o 2 /4 + « 2 )Gq jbi + i8ira 3 hi^3q] /4 -E- ir?G^. b . - 37ra 2 r s - {Un/VS - 16)wa 3 

3 

+ ia 2 (a-r s )^ 3 (27r) 3 (5(q, -b,)]} + ^ £ [ - 8 7 raG q 3 . bi - ia/i J (27r) 3 ( 5(q, - b,)] , (130h) 



where i 7^ i, q- = q + 2b^ (nof the outgoing momenta), y/z 
is defined with a branch cut along the negative z-axis, and 



•^qb 



( 2 + q-b + 6 2 -£-ir 7 )- 1 , 



(131) 



T(- 1 )(b 1 b 2 b 3 ;qiq2q3)= [ - 5127r 3 a 3 cf (q, , b, 



t.j=i 



mn 2 a 3 ( Pl +h t )GZ ]hz ], (135b) 



cf(q,b) = J ^G^GL 



3 ^qk^kb, 



(132) 



cf (q, b, k) = lim 



K^oo 167T 3 

2cf(k,b)- 



w K 
m — 

K 



d 3 k 



k<K 



(27T) 



3 U qk 



^kb 



47T 



(133) 



The loop integrals cf and cf emerge from the Z-S expansions 
(Appendix B) of T^ /2+k _ q/2 _ k and T^ /2+k _ q/2 _ k , 
respectively. 

Remarkably, cf (or any loop integral or lattice sum encoun- 
tered in the present paper) is free from uncontrolled ultraviolet 
divergence; this finiteness follows naturally from the rules of 
the Z-S expansion (as can be easily understood from a similar 
but simpler problem, Example 2 in Appendix B). 

Comparing the above results for 'fqjqj.qa 

iq2qa] 

with Eq. (123), we find the same result for the two-boson T- 
matrix as Eq. (121), and the first 3 terms in the low-energy 
expansion of the three-boson T-matrix: 



T(bib 2 b 3 ;qiq2q3) = ^ T (s) (bib 2 b 3 ; qiq 2 q 3 ), 

s=-2 

(134) 

Where T (s) (bib 2 b 3 ; qiq 2 q 3 ) ~ E s/2 (not excluding 
E s l 2 In" 1 E). 



T<- 2 > (bxbabg; qi q 2 q 3 ) = 64^V £ b( , (135a) 



T (0) (b!b 2 b 3 ; qiq 2 q 3 ) = -6 [D - 24irw(7ir/V3 - 8)a 4 

3 

-36ir 2 a 3 r s ] + 647r 2 a 3 ^ [327r 2 acf (q^, b;, lap 1 ) 



+ (r 8 /2 - a)(p 2 + fc?)G£ bi - a Pj hiG* bi 
+ i8Tra(pj +ft i )cf(qj,b i )]. 



(135c) 



The momenta q^'s satisfy the constraint X^=i = ^> in 
addition to X) i=1 Q« = Z)i=i bj = 0, in Eqs. (135). 

The two lowest order contributions to the 3-boson T-matrix, 
as well as the logarithmic dependence on energy at the third 
order [in a 4 cf (q^, b^, | a. | 1 ) ] , are universal in the sense that 
they depend on the scattering length a only, in agreement with 
Ref. [26]. The leading nonuniversal contributions are 



-6(D - 367rVr a ) + 327r 2 a 3 r s (P? + /l ?) G < 



E 

q 3 b s 



i.j = l 



The momentum-independent term in this expression is pro- 
portional to D — 3Qn 2 a 3 r s , while the leading nonuniversal 
contribution to the BEC energy [in Eq. (118)] is proportional 
to D + I2ir 2 a 3 r s . We conclude that they can not be absorbed 
into a single 3-body contact interaction parameter g 3 , unless 
the two-body effective range r s = 0. This disagrees with 
Ref. [7]. 



C. Comparison with Ref. [7] 

If r s = 0, Eqs. (135) then agree with Ref. [7]. To see this, 
we use the power series for the 2-boson T-matrix to expand 



18 



the last term in Eq. (77) of [7] (note that gi2 in [7] is twice /13 
here), and get all the 4 terms in Eqs. (135) of the present paper 
which contain G E b . . The quantity T^ pl of [7] corresponds to 
the first term on the right hand side of Eq. (135b) above; the 
finite contribution to 7^ 1PI in Eq. (80) of [7] corresponds to 
the last term of Eq. (135c) above. The sum of the divergent 
[7], but dimensionally regularized [7], terms in X, IPI and the 
3-body term — [g 3 (k) + 5g 3 (n)] [7] can be expressed in terms 
of c E defined above: 



(87ra) 4 



cms 



E " E 



cf(k;,k'-,K) - g 3 (K) 



where 



cms 



18io[ln(27r) + 1 -7] + V3[2<5' + 9 In 3 - 18] 



64tt 3 



= 0.055 571 793 27. 



Here 6' = £~ [(V3 + «)~ 2 - (2/3 + n)- 2 ]. The 
momentum-dependent terms in Eqs. (135) thus completely 
agree with Ref. [7] at r s = 0. Further matching the con- 
stant terms, we find the relation between the 3-body parameter 
gs(n) of Ref. [7] and the scattering hypervolume D defined in 
the present paper at r s = 0: 



ffadal -1 ) = 6[D - 24irw(7ir/V3 - 8)a 4 ] + 
= 6{D + 977.736 695 a 4 ). 



(87ra) 4 



-cms 



At r s 7^ 0, the discrepancy between Ref. [7] and our result, 
Eqs. (135), disappears if r s is included in the effective field 
theory. The 2-boson interaction vertex becomes [9] 



,/2\ 



i(— 87ra — 2ira 2 r s k 2 — 2ira 2 r s k 



where k c /2 ± k and k c /2 ± k' are the momenta of the two 
bosons before and after the interaction, respectively, k 7^ k' if 
a virtual particle is involved. The tree diagram contribution to 
the 3-boson T-matrix, Fig. 5. (a) of Ref. [7], is modified as 

3 

E Gqjhi \- 8na + 27Ta2r sP 2 j + 27ra 2 r s (b l + q^/2) 2 ] 



E 



x [87ra + 2na 2 r a h 2 + 2Tro 2 r s (q j + hi/2) 2 

3 3 

- 2887r 2 a 3 r s + 0(E 1 ). 



r s corrections to all other diagrams are 0(E S ), s > 1/2. 
Comparing the modified EFT results with Eqs. (135), we get 



93{m 



6h 2 



(D + I2ir 2 a 3 r s + 977.736 695 a 4 ), 



^boson 

where the SI units have been restored. 



(136) 



D. Implications for the BEC energy and other properties 

The BEC energy per particle at zero temperature was com- 
puted to order p 2 by Braaten and Nieto [7]. The result, 
Eq. (96) of Ref. [7], is expressed in terms of gz{n). 

Although the three-boson scattering amplitude receives r s 
corrections as shown above, the BEC energy per particle does 
not suffer from r s corrections at order p 2 in the thermody- 
namic limit, if it is expressed in terms of #3 (k) [9]. 

Using Eq. (136), we can now express Braaten and Nieto's 
result for the BEC energy [7] in terms of the scattering hyper- 
volume D defined in Eq. (3). The outcome is precisely the 
same as Eq. (118), and C' E in Eq. (1 18) is found to be 

C E = 977.736695/12tt + 8w [ ln(167r) + 0.80 ± 0.005] 



= 118.65 ±0.10. 



(137) 



The present author did an independent calculation of C E , 
using finite-range interactions (details to appear elsewhere), 
and found 

C E = 118.498 920 346 444 (exactly rounded), (138) 

which is nearly the same as Eq. (137). 

For a dilute Bose gas of hard spheres, using Eqs. (4), (118), 
(138) and r s = 2a/3, we obtain the following result for the 
constant defined in Eq. (1): 



C 3 



167.319 69 ±0.000 06. 



(139) 



This completes our calculation of the ground state energy per 
particle of a dilute Bose gas of hard spheres, through order p 2 . 

For a dilute Bose gas with a 3> r e (r e is the range of inter- 
particle forces), Braaten, Hammer, and Mehen found that £3 is 
near 141 [8] and has a small imiginary part associated with 3- 
body recombination [8]. In comparison to this system, a hard 
sphere Bose gas with the same scattering length a, number 
density p, and boson mass, has an energy density that is larger 
by a relative fraction w 26pa 3 , a pressure that is larger by a 
fraction w 52pa 3 , a speed of sound that is faster by a fraction 
w I x 26pa 3 = 39pa 3 , and a specific heat that is smaller by 
a fraction w | x 26pa 3 = 117pa 3 . Here the specific heats are 
compared at the same temperature T <C p 1//3 aT c , where they 
are dominated by phonons with wavelengths 3> 1 / ^fpa. (T c 
is the critical temperature.) 

These differences, ie nonuniversal effects, must extend to fi- 
nite temperatures, including both T <T C and T > T c . (Their 
magnitudes will change when T is raised.) 

For T/T c - O(l) and pa 3 , pr 3 < 1, the thermal de 
Broglie wavelength greatly exceeds r e and a, so the interac- 
tion is well described by the constants a, r s , and D, etc. 

The nonuniversal corrections to T c , as an extension of the 
calculations reviewed in Ref. [27], are of particular interest. 
Is the leading order nonuniversal correction determined by 
r s , or D, or both? The present author speculates it is perhaps 
by r s alone, and is of the form c'pa 2 r s T°, where Tj? is the 
critical temperature of the noninteracting Bose gas, and d is a 
universal numerical constant. 
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VI. SUMMARY 

We have shown that the three-body force near the scatter- 
ing threshold can be predicted in a way very similar to the 
two-body force, ie, by solving the Schrodinger equation at 
zero energy and matching the solution to the asymptotic for- 
mula for the wave function at large relative distances. This 
approach is applicable to the n-body force as well. 

Although in this paper an explicit calculation of the three- 
body force is given for hard spheres only, one can apply the 
general formulas, Eqs. (45), to many other interaction poten- 
tials, to predict the effective three-body forces. The unknown 
wave function has 3 independent variables only, because of 
the translational and rotational symmetry, and is not very dif- 
ficult to study on a present-day computer. 

The author believes that the asymptotics of the three-boson 
wave function at large relative distances found in this paper 
is also applicable to composite bosons with finite-range in- 
teractions. For weakly bound dimers of fermionic atoms, the 
two-dimer scattering length a d w 0.6a/ [28], where a/ is the 
atomic scattering length, but the parameters r s and D for these 
dimers are still unknown, which must be computed to deter- 
mine the equation of state of many such dimers [28-30] more 
accurately [31]. 

We have expanded the ground state energy of three bosons 
in a large periodic volume to the order L~ 7 , using a pertur- 
bation procedure that resembles the derivation of the small- 
momentum structure of <f>^ . The result, Eq. (5), may be 
combined with a Monte-Carlo simulation (for nonrelativistic 
particles such as atoms) or Lattice QCD simulation (for low 
energy pions) to determine the three-body force. 

Equation (5) suggests that, to determine D accurately, L 
should greatly exceed 17a in these simulations, or the higher 
order corrections can overwhelm the lower order terms. 

To extract D from N-body simulations in volumes of mod- 
est sizes, one may find it helpful to derive a systematic expan- 
sion of E at large L, but fixed Ma/ L, corresponding to the 
intermediate regime discussed in Sec. IV E 2. 

For pions, the Compton wavelength is comparble to a and 
r s ; whether the relativistic corrections modify the form of 
Eq. (5) is a question of interest. 

We have computed the scattering amplitude of three bosons 
at low energy to 0(E°), using finite-range interactions. Our 
result, Eqs. (135), disagrees with an EFT prediction [7] at 
r s 7^ 0. By including r s corrections in the EFT, however, 
one can eliminate the discrepancy completely. The resultant 
relation between 53 (k) and D, combined with the many-body 
energy formula of Ref. [7], and the three-body force computed 
in the present paper [Eq. (4)], solves a longstanding problem 
of theoretical interest, namely the complete second-order cor- 
rection to the ground state energy of a dilute Bose gas of hard 
spheres. The result shows small differences between this sys- 
tem and a dilute Bose gas with large scattering length consid- 
ered by Ref. [8]. 

We have studied the energies and momentum distributions 
of Af bosons in a large volume, from which we deduce the 
energies and condensate fractions of dilute BECs with sizes 
greatly exceeding the healing length. The result for the en- 



ergy, Eq. (118), agrees with the effective field theory, but the 
condensate depletion, Eq. (119), disagrees. It is found that at 
low density, even the population of the zero-momentum state, 
or the off-diagonal long-range order [32], is affected by short- 
range physics missed by the EFT. 

Although the three-body force is usually a small correc- 
tion to physical observables in a dilute BEC, dramatic effects 
may be obtained near a three-body resonance [8, 33]. Alter- 
natively, one may reduce the two-body force by tuning the 
scattering length to a zero crossing [34]. At a = the BEC 
ground state energy density « H 2 Dp 3 /6m. 

The author would like to draw the readers' attention to a 
few mathematical techniques. The Z-functions and Z-5 ex- 
pansions (Appendices A and B), an extension of the famil- 
iar 6 function method, facilitate the derivation of the small- 
momentum and the large-volume expansions. The method of 
tail-singularity separation, as described in Appendix C, en- 
ables us to evaluate many lattice sums with virtually arbitrary 
precision. The utility of these methods is certainly not limited 
to the present work. 



Acknowledgments 

The author thanks N. Barnea, J.-R Blaizot, A. Bulgac, W. 
Detmold, Y. Nishida, and M. J. Savage for useful discussions. 
This work was supported by Department of Energy grant DE- 
FG02-00ER41132. 



APPENDIX A: THE Z-FUNCTIONS 

We define generalized functions Z/k 2n and Zb(k)/k 2n+1 
(n= 1,2,3,..., and b > is a constant): 

/ d3k ^P {S) ( k ) = (s < 2n - 4), (A2) 

Jail k K 

I ^ k T^V (2n ^ ) (k) = (fc > and n > 2), (A3) 

J\k\<k k 

I ^iSSi^W - (, < 2n - 3), (A4) 
/ d 3 fcg§^- 2 )(k) = 0. (A5) 

where p( s )(k) is any homogeneous polynomial of k with de- 
gree s (= 0, 1, 2, . . .). Z/k 2 can be identified with 1/fc 2 . 

We can use ordinary functions to approach the Z-functions. 
For instance, (fc 2 - 3r/ 2 )/(fc 2 + rj 2 ) 3 — > Z/k A when 77 — > 0+. 

The Z-functions, like the <5 function, are merely Fourier 
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transforms of some ordinary functions (n = 1, 2, 3, ... ): 

d 3 fc Z ilr .„ (-l)«+l r 2«-3 



(2tt) 3 fc 2i 



_ e ikr = 



47r(2n-2)! 



lim 

r;->0+ 



/ d3 



n+l„2n— 3 



47r(2n-2)! 



-ik-r— r/r 



A: 2 ™ ' 



d 3 fc Z b (fc) c ; k . r = ^(n)^ 



(2tt) 3 fc 2 «+ 4 
lim / d 3 r J F 1 b (n) (r)c- ik r - I ' r = 



(~l) w r 



n„,2n-2 



£2n+l ' 
2n-l 



(A6) 
(A7) 
(A8) 
(A9) 



2?r 2 (2n- 1)! 



ln(6r) + 7 - ^ t] 



i=i 



where 7 = 0.5772 • • • is Euler's constant. 
Zb(k)/k 2n+1 has the following properties: 

Z b (k) Z b ,{k) 47rln(6'/6) T72n _ 2 
^+T-l^Tr= (2n _ 1}! V fe «(k), (A10) 

Z 6 (cfc) _ 1 Z 6/c (fc) 



( c k)2n+l c 2n+l £2n+l 



(constant c> 0). (All) 



APPENDIX B: Z-S EXPANSIONS 

Example 1. Consider the expansion of (|k| 2 + 3<7 2 /4) _1 
at small q. Naively, it is 1/fc 2 - 3q 2 /(4fc 4 ) + 0(q 4 ). But 
this series fails around k = 0. In particular, i? (k, q) = 
(fc 2 + 3q 2 /4) _1 - 1/fc 2 has a finite integral over all k: 

/ ^i? (k,q) = -&q, but the integral of -3g 2 /(4fc 4 ) 
over all k is infinite. 

Now subtract -|^(27r) 3 <5(k) from Ro(k, q) to obtain 
i?i(k,q). Clearly i?i(k,q) w -3<r7(4fc 4 ) at fc > and 
q -> 0, but / i?i(k, q)d 3 fc = 0, so actually 

i?i(k,q)«- — ^, 

where Z/fc 4 is defined in Appendix A. Further subtract- 
ing — 3<7 2 Z/(4fc 4 ) from i?i(k, q), we get a remainder w 
( % /3/647r)g 3 V^(2 7 r) 3 5(k). Continuing this subtraction pro- 
cedure to higher orders in q, we get Eq. (31a). 

The second method to derive the Z-5 expansion of (fc 2 + 
3<; 2 /4) -1 is: first Fourier transform it for fixed q, to obtain 
cxp(— \/3qr/2) j (4nr), then expand it in powers of q: 



_L V3g 3g 2 

4?rr 8?r 327r r 



64tt 



2 , 3 <? 
r H r — 

5127T 



3 ^K r ^ +0{q% 

5120tt W ; ' 



and finally transform the series back to the k-space term by 
term. With the help of Eq. (A7), one gets Eq. (31a). 

Although the above two methods are equally valid, the first 
one is more useful, because the Fourier transforms of most 
functions of k in this paper can not be found analytically. 

Example 2 [needed in deriving Eq. (35)]. At small q 

(|k + q/2|- 1 + |k - q/2|- ^/(fc 2 + 3q 2 /4) 

= 2Z K (fc)/fc 3 + [87r-47r 2 /S/3-87rln(g//s)]5(k) + 0(q 2 ), 



where k > is arbitrary. Higher order corrections may easily 
be obtained as well (see the general rules below). 

Generics. Consider a function F(k, o) (where o are a set 
of variables), with a finite integral over any finite region of 
the k-space for o / 0, and a unique singularity at k = 
for o^0. The Z-S expansion is of the form F(k, o) = 
F z (k, o) + F$(k, o), where F z is a series including ordinary 
and/or Z functions of k and may be directly inferred from the 
Taylor expansion of F(k, o) at small o, and 

F s (k,o) - ^ C „ y (o)V 2m g(^(V fe )(2 7 r) 3 ( 5(k), (Bla) 

mj 

c m3 (o) = (-lJ'^V^Q^CVO^Q^W]}" 1 

x J -^fc 2m Q«(k)[F(k, )-F z (k, )]. (Bib) 

Here are all the independent homogeneous harmonic 
polynomials, satisfying J Q^\k)Q^ )(k)d 2 k = for j ^ 
j'. The degree of is lj. 

If F(k, o) is rotationally invariant around an axis q, only 
are needed above, and the coefficient before the integral 

sign in Eq. (Bib) becomes (2 JjTt^m+ti+i)!! ■ 

If F(k, o) is a completely symmetric and even function of 
k x , k y , and fc z , only those harmonic polynomials with the 
A+ symmetry are needed: Q^(k) = 1, Q {a) (k), Q^(k), 
g (8) (k), where 

Q<»>(v) =vi+v* + v 4 z - 3(v 2 x v 2 y + v 2 y v 2 z + v 2 v 2 x ), (B2a) 
QW( v ) =v 6 x +vl + vl (lh/2)(vtv 2 + V y z + v\v 2 x 



1 24, 24, 2 4\ , fin 2 2 2 

+ V x V y + V y V z + V z V x) + 90v x V y V zi 



(B2b) 



Q 



U(v e x v 2 y +v 6 y v 2 z +v 6 z v 2 x 



+ v>; + K< + «)+ ^(v>l + v y vt + v*X). (B2c) 



APPENDIX C: LATTICE SUMS 

The following methods are used to evaluate lattice sums: 
tail-singularity separation, Poisson summation formula, and 
convergence acceleration (based on the large-n asymptotics 
of the summands). 

To evaluate a lattice sum ^2 n X(n) (sum over 3-vectors 
of integers), where X(n), as a continuous function, has both 
singularity in the real-n space and a power-law tail at large n, 
we sometimes break X(n) in two pieces: X(n) = Xi(n) + 
X2 (n), such that X\ (n) has singularity but no power-law tail 
(ie, decays much more rapidly at large n), while X 2 (n) has 
power-law tail but is sufficiently smooth. E n Xi(n) is done 
directly, while E a ii n -X^n) is approximated by J d 3 71X2(11) 
to a very high precision. We call this method tail-singularity 
separation (TSS). 

For instance, to compute a s (s — 1,2,3- ••)> we write 

n- 2s = E^o ^ CXP ^r 2) + X 2 (n), where V > is 
small. X 2 (n) is very smooth, so Eaiin^2(n) is approxi- 
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mated by an integral. Straightforward algebra yields 



a s = c s ir ' T] 



3/2 s-3/2 _ V_ 



Si 



+ EE 

n^O i=0 



— 1 (r/n 2 ) 4 e-"™ 2 



i! n 2s 



(CI) 



where ci = -2, c 2 = 2, c 3 = 1/3, c 4 = 1/15, c 5 = 1/84, 
eg = 1/540, and the error <~ 0(e~ 7r results from 
the approximation ^ all n X 2 (n) w J" d 3 nX 2 (n) [35]. Equa- 
tion (CI) applies to both s = 1 and s = 2, 3, • • • . At 

= 1/10, one already gets about 43-digit precision for a s . 

The TSS is not an arbitrary exponential accelaration 
method. For example, J2 a ^o 71 3 * s not sufficiently acceler- 
ated by 2^ n ^o exp(— ?7n 3 )n _3 , because [1 — exp(— rjn 3 )]n~ 3 
is not smooth enough: it contains a term <~ |n| 3 at small n 
which is singular. Inspired by Takahasi and Mori's quadrature 
method [36], we obtain a proper TSS formula: 

1/n 3 - nrj/2n 2 = 2/[l + exp(7r sinhr/n)]n 3 + 5(n), 

where S(n) is smooth. We then easily derive (all the lattice 
sums below, except /3 1A , are used in computing Co and C\): 



a 15 = ji + 47rln?7 + irair]/2 + tt(tt 2 - 2)ry 3 /24 
+ ^ 2/[l + cxp(7r sinhryn)]™ 3 + 0(e-* 2/r >), (C2a) 

f R /7T 

ji= lim 47r / x^ 1 tanh ( — sinha; )dx — AtrlnR 
R^oo J \2 J 

= 16.489 380 548 112 915 838 168 866 783 965 159 811. 
Similarly 

a 2 .5 = Ka 2 'q/2 + j 2 f] 2 - tt(tt 2 - 2)a 1 T] 3 /24 
- 7r(7T 4 - 5tt 2 + 1)^/240 + ^ 2/[l + cxp(7rsinhr/n)]n 5 

n^O 

+ 0(0-^/"), 



(C2b) 



/■oo 

j2 = 47T / x~ 3 tanh ^— sinhxj — ir/2x z 
= -24.436 776 868 803 521 072 197 485 676 562 043 398. 
Let 

ai 9) = lim V c-" n Q (9) (n)/n 2s , (C3) 



n/0 



and similarly for a« and ai 8 . Here Q^, (J' 4 ' and are 
defined in Eqs. (B2). 

For any integer s > 1, we derive a TSS formula 

"i^E^S^-'+Ote-""). (04) 

and similarly for and ai 8 - 1 . 



Using the Poisson summation formula and the TSS method, 
we get 

a { 5 9 l = ~(2\ 5 /9l\)Y / Q {9 H^)Kvn) + 0(e^ 2 ^), (C5a) 

n^0 

a« = (2V/13!!) y Q('»(n)(W + 0(e~^^), (C5b) 

4 8 1 = -(2 10 7r 9 /17!!) y Q< 8 )(n)iW+0(o- 2 /'), 

(C5c) 

where l(x) = In [tanh sinhx)] . 

Using the Poisson summation formula, we get 



W A (n) = -V3 7T 2 n + 7T y e _V5 



nmn— lirm-n 



m^0 



llfl(n) = (27r 2 /V3n) J] c-^ 3 



7rmn— i7rm-n 



/m, (C6a) 



(C6b) 



We derive TSS formulas 



6>Ai(n) = -r?(l - c-^ n2 )/n 2 + 2ir dm 
x (1 _ e -^ Bmc)(1 _ e -^ 2 )/5„ mc + y 



m 



m/0 



x (e^ 5 — + e - T ' m2 -e-'' s '™- T '' n2 ) + 0{c-^ l 2r] ), (C6c) 
6» B i(n) ^/[(^^n 2 )^''" 2 - l]/n 4 + 27r / dm dc 

x [1 - (1 + V S nmc )e-^} (1 - e-"™ 2 )/5 2 mc 
+ E S nV 2 [(l+lMe^" m (l-e^ m V^ ,m2 ] 

rn^O 

+ 0(c-^ 2/2 ^), (C6d) 

where S nm = n 2 + n- m + m 2 , and SVimc = n 2 |nmc|m 2 . 
With Eqs. (C6), p A1 (n) and p B1 (n) [Eqs. (86a) and (94a)] are 
evaluated very accurately for any finite 3-vector of integers 
n ^ 0. 

In addition, at large n we derive asymptotic formulas 

p A1 (n) = p^l ) (n) + 0(n- 12 ), (C7a) 
/^(n)^ ir 2 w/n+ 2ai/n 2 + 4/3n 4 +4a^ ) Q (9) (n)/15n 10 

+ 4af ) Q«(n)/231n 14 + 4a^ 8) QW(n)/195n 18 , 
p sl (n)=pg 1 3) (n) + 0(n- 14 ), (C7b) 
p£?(n) = 2V3 7r 2 /n 3 +2 ai /n 4 + 2/n 6 +4a< 9) Q(^(n)/3n 12 

+ 4a^g w (n)/33n 16 + 12a^ 8) QW(n)/65n 20 . 
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For a 1A1 , a 2A1 , and a 1B1 , we thus have 

a sAi = K 2 wa Q5+s + 2 ai a 1+s + 4a 2+s /3 + 4a^ ) a^ s /15 
+ 40^4^/231 + 44 8) 4t/195 
+ E n_2 *[^i(n)-Pi 1 i 1) (n)], (C8a) 

a lm = 2V3n 2 a 25 + 2a l( x 3 + 2a A + 4a[ g) af ] /3 
+ 4a^ ) 4 i) /33 + 124 8) 4i } /65 
+ E^ 2 [^i(n)-pi, 1 1 3) (n)], (C8b) 

n^O 

and the sums over n are truncated at \n x \, \n y \, \n z \ < 15 to 
yield results for a 1A1 , a 2A1 , and a 1B1 with more than 18- 
digit precision. 

From Eq. (99) we deduce 

&1AA1 — ~ a \AA\ ~~ 7i" 3 w(77r/\/3 — 8)a 1 + 8ir 3 wa 1 ln(27r) 

— a 1 a 1A1 + 3a 1A2 + 3a 1B1 — {a\ + a 2 )a 2 + 6a x a 3 — 9a 4 , 

(C9) 

a iAAi = J im E n ~ 2 PAAi( n ) - STT 3 w ai InN 

N—>oo z — ' 

n^O; n<N 

+ 47r 4 wiV(81n7V- 16 + 77T/V3). (CIO) 



It can be shown that a 1A2 = a 2A1 , and 



a 1AA1 = lim p A1 (n) - Att 5 w 2 N - IQ^wa^ In N. 

N^oo — ' 

n^O;n<AT 

(Cll) 

At large n we derive from Eq. (C7a) 



p^(n)=^ 12 )(n) + 0(n- 13 ), 



(C12) 



K {12) {n) = tt 4 w 2 /ti 2 + A^wax/n 3 + Aa\jn A + 8n 2 w/3n 5 + 16ai/3n 6 + 8n 2 wa[ 9) Q i9) (n)/15n n 
+ [l6/9n 8 + 16aia^Q (9) (n)/15n 12 ] + 8ir 2 wa { * ] Q^ (n)/231n 15 + [32a^Q^(n)/45n 14 + 16aia^Q w (n)/231n 16 ] 
+ 8w 2 wa[ 8 (n)/195n 19 + [64a^ )2 /4725n 12 + Ma[ 9)2 Q^ (n)/3575n 16 + (l28a^ )2 /10395 + 32a? } /693)Q W (n) /n 18 

+ (l6a< 9)2 /585 + 16a 1 af ) /195)g (8) (n)/n 20 ]. (C13) 

So 

a 1AA1 — n 4: w 2 a 1 + A / n 2 wa l a l b + Aa\a 2 + 8 / n 2 wa 2 5 /3 + 16aia 3 /3 + 8ir 2 wa^ + 16a 4 /9 + 160^4^4^/15 
+ 87r 2 wa? ) a^ ) 5 /231 + 32a[ 9) alf ) /A5 + IQonaf af /2?>l + 8ir 2 wa {8) a { gl/ 195 + 64a^ )2 a 6 /4725 + 64a^ )2 4 9) /3575 
+ (l28a^ )2 /10395 + 32a? ) /693)4 ) + (l6a< 9)2 /585 + 16ai4 8) /195)4o + ^ [p\M) - K {12) {n)] , (C14) 

I 



where the sum over n converges very rapidly; truncating it 
at \n x \, \n y \, \n z \ < 15 we obtain the result for a 1AA1 with 
18-digit precision. 



1-fold lattice sums (rounded to fit the 2-column format): 



The results for the above lattice sums are listed below. They 
determine Co and C\. 



ai = -8.913 632 917 585 151 272 687 120 136, (C15a) 

ai. 5 = 3.821 923 503 940 635 799 730 123 034, (C15b) 

a 2 = 16.532 315 959 761 669 643 892 704 593, (C15c) 

a 2 . 5 = 10.377 524 830 847 083 864 728 948 355, (C15d) 
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a 3 = 8.401 923 974 827 539 993 146 138 987, (C15e) 

a 4 = 6.945 807 927 226 369 624 170 778 023, (C15f) 

a 6 = 6.202 149 045 047 518 551 930 416 392, (C15g) 

a[ 9) = 1.127 757 148 686 792 075 014 583 731, (C15h) 

a[ 9 \ = 5.672 605 625 422 259 129 524 572 370, (C15i) 

4 9) = 5.772 772 158 341 296 181 095 291 055, (C15j) 

c4 9) = 5.890 866 300 404 517 457 312 864 285, (C15k) 

a { 8 9) = 5.947 443 615 615 246 912 541 646 623, (C151) 

= -2.434 385 049 385 522 287 574 679 853, (C15m) 

a W = 5.242 702 841 443 828 214 862 689 624, (C15n) 

af = 5.451 072 910 162 494 576 739 294 958, (C15o) 

af> = 5.715 464 651 152 249 495 454 861 607, (C15p) 

af ] = 16.571 410 717 493 131 178 469 668 510, (C15q) 

afl = 6.156 579 477 795 506 959 879 789 010, (C15r) 

a$ = 6.109 685 927 989 065 132 252 619 095, (C15s) 

a$ = 6.054 088 533 916 337 242 283 438 742. (C15t) 



2- fold lattice sums: 

a 1A1 = -190.172 897 984 865 754 8, (CI 6a) 

= = 111-807 832 628 721 133 609, (C16b) 

a 1B1 = 221.523 005 657 695 107 22, (C16c) 

3 - fold lattice sums: 

a 1AA1 = -2996.889 395 378 764 86, (C17a) 
&1AA1 = -6591.229 842 103 439 89. (C17b) 

Finally, from Eq. (C6a) we get 

o , /q 2 „ i ^ \ ^ „— \Z3irmn— i7rm-n /„2 rt _ 1 

p\A = — vo7T a 05 +7r > e /n m 

= 48.614 754 175 227 821 038 934 419 912. (C18) 

(One can show that a Q 5 = a 1 /ir.) (3 1A is needed in 
Eq. (HOf). 
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